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Previously  reported  studies  of  the  HgCdTe/Photox™  SiC^  interface  on  MIS 
devices  have  shown  that  the  net  charge  of  this  interface  depends  on  pre-passi¬ 
vation  surface  treatments,  water  content  and  surface  fermi  energy.  Studies  in 
this  report  have  predicted  that  the  charge  contribution  due  to  interface  traps 
will  be  different  for  n  and  p  type  material  and  demonstrated  this  effect  on 
both  diodes  and  MIS  devices.  Admittance  spectroscopy  and  quasistatic  C-V 

measurements  have  shown  the  slow  state  contribution  to  charge  to  be  due  to 

_2 

states  communicating  with  the  oxide  in  the  10  Hz  range.  Charge  instability 
of  the  interface  was  demonstrated  to  be  due  to  charge  injection  in  contrast  to 
mobile  ions  as  in  the  Si/Si02  system.  Capacitance-voltage  measurements  at  low 
temperature  have  revealed  a  discreteness  of  inversion  layer  formation  or  trap 
levels.  Pre-passivation  surface  treatments  have  a  serious  impact  on  interface 
traps  as  well  as  bulk  trap  characteristics.  Polished  and  etched  surfaces  have 
shown  a  redistribution  of  a  mid-gap  bulk  trap  level  as  seen  by  deep  level 
transient  spectroscopy.  The  c  ncentration  of  this  trap  will  strongly  affect 
G-R  current  and  diode  RoA. 

The  implantation  task  has  shown  that  electrical  activity  due  to  implant 
damage  follows  the  random  portion  of  the  atomic  profile,  but  raised  an  order 
of  magnitude,  for  both  random  and  channeled  samples  of  B,  P  and  Cl.  Channeled 
hydrogen  shows  distinctly  n  type  damage  but  with  no  direct  correlation  to  the 
channeled  profile  obtained  from  SIMS.  Damage  in  randomly  implanted  samples 
was  found  to  create  a  weakly  repulsive  trap  at  =>Ec-0.20  eV  (Eg=0.27eV)  with  a 
large  barrier  height.  The  concentration  of  this  trap  scales  linearly  with  the 
implant  species  concentration.  Continuing  use  of  angle  resolved  photoemission 
spectroscopy  to  observe  the  surface  and  bulk  electronic  band  structure  has 
been  applied  to  samples  with  X=i.O,  0.31,  0.23  and  0.00.  The  dispersion  and 
polarization  dependences  seen  are  similar  to  those  for  CdTe.  There  is  some 
divergence  of  theory  and  data  but  the  general  agreement  validates  the 
assumption  of  a  nearly  free  electron  final  state. 

Calculations  showing  unusual  bond  length  relaxation  characteristics  of 
HgCdTe  predict  a  weak  HgTe  bond  which  leads  to  the  mechanical  instability  of 
the  lattice.  This  simple  radial  force  model  agrees  with  the  trend  predicted 
by  Mikkelson  and  Boyce  for  InAs-GaAs.  The  effects  of  this  unusual  bonding 
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behavior  on  the  thermodynamics  and  kinetics  of  interfaces  has  been  studied. 

Ion  milling  of  HgCdTe  LPE  heterojunctions  resulted  in  HgCdTe  exhibiting  a  more 
metallic-type  recovery  than  CdTe.  This  results  in  lower  defect  densities  and 
narrower  stacking  faults.  This  also  supports  the  theoretical  predictions  of  a 
weakened  HgTe  bond  and  ease  of  Hg  vacancy  formation. 

Soft  x-ray  photoemission  spectroscopy  was  used  to  monitor  the  kinetics  of 
activated  oxygen  uptake  on  the  cleaved  (110)  face.  The  oxides  grown  were 
composed  of  Cd  and  Te  oxides.  Hg  was  shown  to  be  lost  from  the  surface  during 
the  initial  stages  of  oxide  formation  but  no  Hg  was  found  in  the  final  ox¬ 
ide.  Monitoring  of  the  growth  rate  on  samples  of  various  x-values  showed  the 
growth  of  the  first  two  monolayers  to  scale  with  the  number  of  Hg-Te  bonds. 

Using  an  extension  of  Harrison's  method,  dislocation  energies  per  unit 

-3  -9 

length  were  calculated  to  be  proportional  to  d  to  d  and  hardness  propor¬ 
tional  to  d-^  to  d"*1  where  d  is  the  bond  length.  These  calculations  are 
confirmed  by  experiment.  Hardness  was  found  to  be  related  to  the  interaction 
energies  among  dislocations.  It  is  proposed  that  an  alloy  with  a  second 
constituent  of  a  shorter  bond  length  will  improve  crystal  strength.  Experi¬ 
mental  verification  of  this  has  been  seen  on  GaP  in  GaAs  and  4%  Zn  in  CdTe. 
Calculations  show  the  already  weak  HgTe  bond  to  be  further  destabilized  by 
adjacent  CdTe,  HgS  or  HgSe  but  to  be  stabilized  by  ZnTe.  This  leads  to  the 
conclusion  that  HgZnTe  will  have  a  structural  advantage  over  HgCdTe. 
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This  report  summarizes  results  obtained  on  the  study  of  the  HgCdTe  sur¬ 
face  for  the  six-month  period  ending  July  1,  1984.  The  purpose  of  the  study 
is  to  develop  a  microscopic  picture  of  the  surface  and  near  surface  regions  of 
HgCdTe,  particularly  the  defect  structure,  which  play  a  direct  role  in  the 
surface  passivation  of  IR  detectors. 

Successful  passivation  of  detectors,  particularly  photovoltaics,  is  of 
critical  importance  for  the  smaller  size  elements  found  in  large,  high-per¬ 
formance  arrays.  Here  the  volume  to  surface  area  ratio  has  decreased  suffi¬ 
ciently  so  that  it  is  no  longer  adequate  to  optimize  bulk,  junction  parameters 
to  achieve  high  performance.  Phenomena  determined  by  the  surface  properties 
come  to  completely  limit  device  operation.  These  effects  are  not  predictable 
from  bulk  diode  modeling  and  it  is  now  necessary  to  design  and  model  surface 
properties  in  conjunction  with  bulk  properties.  These  surface  properties 
cannot  be  inferred  from  bulk  properties  but  must  be  measured  separately  in  a 
system  which  incorporates  the  passivation  to  be  used. 

The  effectiveness  of  HgCdTe  passivation  with  Photox”  Si02  has  been  shown 
to  depend  most  strongly  on  the  interface  charge.*  Under  the  right  conditions, 
this  charge  can  induce  significant  leakage  current  at  the  surface  in  a  region 
that  would  otherwise  be  depleted  of  carriers,  effectively  shunting  the  bulk 
junction  RQA.  The  effects  of  interface  charge  on  the  surface  of  an  n+  on 
p  junction  is  schematically  represented  in  Figure  1.  The  transition  from  n  to 
p  type  and  hence  the  junction  is  represented  by  the  vertical  solid  lines  with 
the  depletion  region  limits  for  reverse  bias  shown  by  the  dashed  lines.  The 
p  side  depletion  region  is  wider  since  this  is  the  lower  doped  side.  This 
figure  illustrates  in  (a),  that  for  net  positive  interface  charge  the  p  side 
is  depleted  at  the  surface  while  the  n  side  is  accumulated;  and  (b),  for  net 
negative  charge  the  p  side  is  accumulated  while  the  n  side  is  now  depleted. 

For  the  large  difference  in  doping  from  n+  to  p  sides  as  is  common  for  im¬ 
planted  junctions,  for  example,  the  depletion  region  on  the  n+  side  is  very 
much  narrower  so  that  when  this  side  is  accumulated  the  effective  decrease  in 
the  surface  junction  depletion  width  is  insignificant  in  comparison  to  accumu¬ 
lating  the  p  side.  Pinching  off  the  depletion  region  at  the  surface  for 
negative  net  interface  charge,  in  this  case,  can  increase  the  surface  trap 
tunneling,  in  the  presence  of  appropriate  surface  states.  This  simplest  and 
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Schematic  illustration  of  the  effects  of  interface  charge  on 
surface  of  n+  p  junction. 


straightforward  scenario  does  appear  to  be  a  significant  mechanism  of  surface 
leakage  in  many  MWIR  double  layer  heteroj unction  (DL.HJ)  devices.  There  are 
however  many  complications  on  this  once  we  consider  implanted  junctions,  LWIR 
devices  or  even  the  more  complicated  barrier  junctions  which  are  produced. 

Since  net  interface  charge  is  contributed  to  by  fixed  insulator  charge 
and  charged  interface  states  in  about  equal  proportion,  it  is  not  at  all  accu¬ 
rate  to  model  surface  leakage  based  on  the  same  charge  affecting  both  the 
n  and  p  sides.  The  surface  Fermi  energy,  Efs»  is  different  at  n  and  p  inter¬ 
faces  so  that  interface  traps  are  filled  to  different  levels  resulting  in  dif¬ 
ferent  interface  trap  charge  components.  The  actual  value  of  E^s  is,  of 
course,  affected  by  the  field  induced  junction  and  therefore  must  be  deter¬ 
mined  in  a  self-consistent  manner  accounting  for  the  value  of  ,  rather  than 
being  assumed  from  the  bulk  value  of  E^.  In  general  the  interface  trap  con¬ 
tributions  are  expected  to  be  negative  on  the  n  side  and  positive  on  the 
p  side. 
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The  passivated  surface  region  can  be  further  complicated  by  surface  pro¬ 
cesses  which  may  alter  the  composition  at  the  surface  to  a  different  x-value 
or  produce  a  non-stoichiometric  region  by  depletion  of  Hg  and/or  Cd.  This  of 
course  will  affect  atomic  or  crystal  bonding  (and  therefore  interface  state 
density)  as  well  as  doping  density  and  band  gap.  These  latter  also  have  a 
direct  effect  on  Efg.  Surface  damage  due  to  processing  is  of  fundamental  irn- 

9 

portance  to  HgCdTe  because  of  the  low  binding  energy  of  Hg  in  the  alloy. 

This  has  motivated  the  fundamental  studies  of  alloy  hardness  and  bonding  dis¬ 
cussed  in  Section  IV. 

It  is  also  of  course  necessary  that  the  passivation  layer  be  a  good  insu¬ 
lator;  have  good  adhesion  and  be  chemically  inert.  These  latter  properties 
are  not  considered  in  this  study  since,  although  practically  important,  they 
are  not  the  principal  determinant  of  surface  electronic  structure. 

To  characterize  the  structure  of  the  interface  region  with  respect  to 
charge  distribution,  electronic  states,  and  crystal  defects,  the  study  (and 
this  report)  are  divided  into  four  principal  tasks:  (1)  electronic  structure 
of  the  passivated  interface  (Photox"  SiC^  on  HgCdTe);  (2)  surface  and  bulk  de¬ 
fects  at  interfaces  with  ion  implanted  HgCdTe;  (3)  electronic  structure  of  the 
bare  surface  (cleaved  in  vacuum)  as  well  as  its  oxidation  chemistry;  and  (4) 
theoretical  modeling  of  the  defect  structure  of  HgCdTe,  particularly  the  Hg 
vacancy,  its  activation  energy,  and  energy  levels. 
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Section  I 
Interface  Study 

This  section  presents  results  of  the  continuing  characterization  of  the 
passivated  HgCdTe  interface.  The  effort  is  focused  on  studying  the  interface 
charge,  interface  trap  structure,  and  alteration  of  the  bulk  structure  at  the 
surface  by  pre-passivation  processing. 

The  net  charge  affecting  the  state  of  the  surface  has  three  contributors 
which  are  each  distinctly  important  due  to  the  difference  in  their  response 
times  to  changes  in  applied  gate  bias.  These  are:  (1)  fixed  charge  (Qfix)  in 
the  insulator  which  does  not  respond  to  bias  changes;  (2)  slow  states  (Qs^ow) 
which  are  viewed  as  oxide  traps  communicating  with  the  interface  and  which  re¬ 
spond  on  the  order  of  1  Hz  or  less;  and  (3)  interface  trap  charge  which  re¬ 
sponds  very  quickly.  These  charges  sum  to  give  the  net  interface  charge  which 
determines  whether  the  passivated  surface  is  accumulated,  depleted  or  in¬ 
verted. 

The  "fixed"  charge  contribution  is  due  to  the  work  function  difference 
between  the  applied  gate  metal  and  the  HgCdTe  substrate  in  addition  to  immo¬ 
bile  ions  in  the  oxide.  These  ions  arise  either  during  formation  of  the  Pho- 
tox™  layer  or  by  diffusion  after  the  layer  is  deposited.  The  Photox™  Si02 
used  as  a  passivation  on  HgCdTe  is  inherently  porous  and  permits  diffusion  of 
ions  in  to  the  oxide.  Fast  diffusing  metals,  such  as  indium,  have  been  seen 
to  alter  the  amount  of  fixed  charge  as  well  as  the  effective  dielectric  con¬ 
stant  of  the  layer.  Titanium  is  commonly  used  for  gates  in  our  devices  be¬ 
cause  it  does  not  show  such  diffusion,  and  hence  does  not  contribute  to  fixed 
charge. 

Net  charge  of  the  HgCdTe/Si02  system  has  been  found  to  be  localized  in  a 
layer  at  the  interface.*  Water,  arising  from  the  deposition  or  from  later  ab¬ 
sorption  of  atmospheric  water,  has  been  shown  to  have  a  large  impact  on  this 
charge."*  Water  is  a  natural  by-product  of  the  silane-nitrous  oxide  reaction 
used  to  obtain  Si 02*  The  amount  of  water  (and  SiOH  and  SiH)  created  will  de¬ 
pend  on  deposition  parameters  such  as  temperature,  gas  flow  rates,  and  gas 
mixture  ratios.  In  a  study  to  determine  the  location  of  charge  in  the  system 
it  was  found  that  less  charge  was  present  in  samples  which  had  15(JU  A  of  Si02 
deposited  in  50UA  increments,  as  opposed  to  a  single  deposition.  In  this  case 
the  incrementally  deposited  sample  sees  a  fifteen  minute  1UU°C  bake  between 
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layers  which  may  minimize  water  content.  It  is  also  expected  that  the  "inter¬ 
face"  between  SiC^  layers  will  act  as  a  barrier  preventing  diffusion  of  water 
to  the  interface.  In  the  case  of  absorption  of  water  from  the  atmosphere, 
possible  surface  states  are  created  at  0.06  eV  and  0.14  eV  which  have  been 
seen  to  contribute  up  to  O.ObV  to  flat  band  bias  measure  by  C-V.  This  water 
does  not  appear  to  make  any  fixed  charge  contribution  since  the  surface  poten¬ 
tial  shift  occurs  abruptly  as  the  surface  states  are  seen.  Water  affects  the 
interface  charge  through  the  trap  charge  component. 

At  interfaces  where  Dit  is  large  enough  that  trap  charge  is  a  significant 
contributor  to  net  interface  charge  it  is  expected  that  the  interface  charge 
will  be  different  at  similar  interfaces  with  either  n  or  p  type  HgCdTe.  The 
density  of  interface  traps,  seen  as  a  U-shaped  distribution  in  the  semiconduc¬ 
tor  band  gap  with  C-V  techniques,  is  believed  to  be  made  up  primarily  of 
scattering  induced  tailing  of  valence  and  conduction  band  states  into  the  band 
gap.  This  is  illustrated  in  Figure  2  where  an  experimentally  determined  dis¬ 
tribution  (solid  line)  is  shown  with  the  band  tailing  contributions  necessary 
to  produce  the  distribution.  States  arising  from  conduction  band  tailing  are 
most  likely  acceptor  like  and  those  from  the  valence  band  donor  like.  (This 
was  incorrectly  reversed  in  the  previous  interim  report).  Recall  that  accep¬ 
tor  like  states  are  neutral  when  empty  and  negative  when  filled  while  donor 
like  over  are  neutral  when  filled  and  positive  when  empty.  There  states  are 
filled  according  to  Fermi-dirac  statistics  so  that  the  expected  total  trap 
charge,  from  those  two  types  of  states  in  a  symmetrical  distribution,  is  nega¬ 
tive  at  an  n  type  interface  and  positive  at  a  p  type  interface.  (This  result 
was  correctly  stated  in  the  previous  report).  The  drop-off  of  D^t  in  Figure  3 
near  each  band  edge  is  not  real  but  an  artifact  of  the  measurement  technique. 
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Wafers  of  HgQ^CdQ^Te  were  cut  from  ingots  grown  by  the  horizontal  zone 
melt  technique.  These  wafers  are  slightly  polycrystalline  and  of  random 
orientation.  Isothermal  annealing  in  Hg  vapor  at  250°C  gave  n-type  material 
with  Np  -  Na  ~  7  x  10^  cm-^.  Capacitors  were  then  fabricated  by  deposition 
of  1500 A  of  PHOTOX"  Si02  and  thermal  evaporation  of  400  A  of  Ti,  followed  by 
4000A  of  Au.  The  metals  were  delineated  by  use  of  a  shadow  mask  during  evapo¬ 
ration.  An  indium  contact  soldered  to  the  HgCdTe  provided  the  ground. 

High  frequency  (1  MHz)  capacitance-voltage  measurements  at  80K  (the  mid¬ 
range  of  temperature  scanned  by  DLTS)  provided  information  on  flatband  volt¬ 
age,  depletion  width  and  accumulation  voltage  needed  to  select  DLTS  measure¬ 
ment  biases.  The  values  measured  were:  flatband  voltage  of  +0.1V  ±  0.05V, 

Nd  -  Na  =  7  x  101*  cnT^  ±  3  x  I01*  cm-^ .depending  on  the  individual  device, 
and  a  maximum  depletion  width  of  ~1 .0  un.  Thus  +1.0  or  +0.5V  was  sufficient 
to  accumulate  the  surface  and  fill  traps  and  -3  volts  or  more  would  create  the 
maximum  depletion  region  width. 

DLTS  was  then  used  to  probe  trap  activation  energies  and  response 
times.  Measurements  were  taken  with  a  signal  frequency  of  10  MHz,  50  mV 
amplitude  and  fill  pulse  duration  of  1  ms.  Peaks  were  observed  over  a  range 
at  temperatures  between  1 5K  and  200K.  Fill  pulse  and  measurement  pulse 
heights  were  varied  between  +1 .0V  and  -4.0V  in  0.25  to  1 .0V  increments  to 
profile  the  depletion  region  (maximum  width  =*1 .0  un) . 

The  trap  emission  rate  (e)  is  given  by: 

e  =  1/t  where  x  =  (t^  -  t7)/ln(t^  -  t0) 
in  which  tj  and  t2  are  the  sampling  times  measured  from  the  end  of  the  fill 
pulse.  Arrhenius  plots  of  e  versus  1000/T  yield  then  the  activation  energy, 

EDLTS* 
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Electron  and  hole  traps  in  HgCdTe  have  been  of  interest  due  to  their 
effect  on  generation-recombination  rates  and  the  attendant  adverse  impact  on 
lifetime  (and,  therefore,  device  performance).  Work  has  been  done  to  charac¬ 
terize  these  traps  for  p-type  Hgj_xCd_xTe  for  various  x-values.*  Two  hole 

traps  were  seen  which  possibly  correlate  with  the  mercury  vacancy  levels. 

2 

Electron  capture  characteristics  of  this  trap  have  also  been  evaluated,  as 

well  as  theoretical  and  empirical  investigation  of  this  trap’s  effect  on  diode 
3 

performance. 

To  date,  work  in  HgCdTe  has  focused  on  trap  levels  in  the  bulk  of  the 
material.  However,  it  has  been  seen  that  defect  formation  at  the  surface, 
particularly  compositionally  significant  cation  removal,  is  unavoidably  asso¬ 
ciated  with  the  chemical  and  mechanical  treatments  used  in  surface  prepara¬ 
tion.^  This  can  cause  a  significant  variation  in  electronic  properties,  in¬ 
cluding  trap  activation  energies,  capture  cross  sections,  and  populations  be¬ 
tween  the  material  surface  and  the  bulk.  There  is  also  the  possibility  that 
shallow  level  traps  in  the  bulk  may  become  deeper  states  at  the  interface  (as 
proposed  by  Allen,  Buisson,  and  Dow"’)  thus  creating  surface  states  with  a 
unique  activation  energy  but  due  to  the  same  impurity.  It  is  important,  then, 
from  the  point  of  view  of  device  surface  passivation,  to  examine  traps  near 
the  material  surface.  Surface  state  populations  and  activation  energies  have 
been  characterized  on  the  Si/SiC^  interface^'^  with  constant  capacitance  DLTS; 
however,  this  has  not  previously  been  done  on  the  on  HgCdTe/Sit^  interface. 

This  work  investigates  for  the  first  time  electron  traps  near  the 
HgCdTe/Si02  interface  using  DLTS  on  n-type  Hg^  7 CdQ  3TB  MIS  structures  in  the 
temperature  range  from  15K  to  200K.  Three  states  were  found  and  character¬ 
ized,  providing  activation  energy,  capture  cross  section,  and  depth  profile 
inf ormat ion . 
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DEEP  ELECTRON  TRAPS  NEAR  THE  PASSIVATED 


INTERFACE  OF  HgCdTe* 


V.  A.  Cotton,  J.  A.  Wilson,  and  C.  E.  Jones 


Santa  Barbara  Research  Center 
Goleta,  CA  93117 


ABSTRACT 

Electron  traps  in  bulk  n-type  HgQ.7CdQ.3Te  were  investigated  near  the 
HgCdTe/SiOj  interface  by  use  of  deep  level  transient  spectroscopy  (DLTS)  on 
Metal-Insulator-Semiconductor  (MIS)  structures [D.V.  Lang,  J.  Appl.  Phys.,  3022 
(1974)].  Three  electron  traps  are  found  with  activation  energies  (relative  to 
the  conduction  band  edge)  of  0.12,  0.172  and  0.079  eV,  and  corresponding  cap¬ 
ture  cross  sections  of  3.5  *  10-^,  1.1  x  10-^,  and  1.2  *  10  ^  cm^.  Depth 
profiles  from  the  surface  to  ~l  .0  un  show  the  concentration  of  the  0.172  eV 
trap  to  be  uniform  while  the  0.12  eV  trap  shows  a  strong  depth  dependence.  It 
is  undetectable  at  the  surface,  rising  in  concentration  to  twice  its  bulk 
value  at  ~0.5  ym  depth  then  falling  to  a  value  comparable  with  that  of  the 
0.172  eV  trap  in  the  bulk  (1  un).  Divalent  trapping  behavior  has  also  been 
detected.  In  this  case,  the  deeper  state  significantly  depopulates,  enabling 
a  shallower  state  to  then  depopulate  [J.S.  Blakemore,  Semiconductor  Statistics 
156  (1962),  Permagon  Press,  London]. 

*  Supported  by  DARPA  contract  MDA-903-83-C-0108. 
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"ideal”  quasistatic  curve  (solid  line)  to  obtain  the  Dit  distribution  of 
Figure  12.  The  inversion  side  of  the  measured  quasistatic  is  widened  compared 
to  the  ideal  resulting  in  higher  values  from  -0.10V  to  -0.17V  of  surface 
potential.  It  is  not  possible  from  this  measurement  to  discriminate  between 
real  interface  states  and  charge  nonuni formi ties .  Admittance  spectroscopy 
(discussed  earlier)  and  deep  level  transient  spectroscopy  (DLTS)  are  not 
usually  affected  by  these  nonuniformities  due  to  their  observation  of  a  speci¬ 
fic  response  time.  When  nonuniformities  are  a  factor  in  admittance  measure¬ 
ments  it  is  easily  discernable  from  other  effects. 

DLTS  is  additionally  valuable  as  a  tool  for  fully  characterizing  trapping 
sites  (activation  energy  capture,  cross  section,  and  trap  concentration)  and 
profiling  their  distribution  in  energy  and  concentration  with  depth  into  the 
material.  This  has  also  been  used  to  evaluate  how  bulk  properties  are  altered 
in  the  region  of  a  processed  surface  or  interface.  Results  of  this  prelimi¬ 
nary  DLTS  study  are  discussed  in  Appendix  I  to  this  Section. 
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only  in  ramping  from  accumulation  to  inversion,  however,  and  a  very  interest¬ 
ing  effect  is  seen  as  we  monitor  the  opposing  ramp  with  temperature.  As  we 
begin  our  observation  at  7K,  there  is  a  wide  depletion  dip  within  which  there 
are  two  very  distinct  peaks.  As  we  increase  temperature  these  peaks  become 
less  pronounced,  the  depletion  dip  is  narrowed  and  deepened  and  by  the  time  we 
reach  90  to  100K  we  have  only  the  familiar  single  dip  representing  a  well 
behaved  quasistatic  curve.  Among  the  possible  causes  for  this  phenomena  are 
discrete  inversion  layer  energy  levels.  More  detailed  analysis  is  underway. 

An  important  consideration  in  all  of  the  C-V  measurements  discussed  is 
the  effect  of  charge  nonuniformities  on  the  device  parameters  thus  extracted, 
most  notably  on  effective  interface  trap  density  (D^t).  When  these  nonuni¬ 
formities  occur  with  distribution  lengths  on  the  order  of  a  depletion  region 

width  (*lum)  of  the  MIS  device  measured,  which  in  our  case  has  an  area  of 
-3  2 

2  x  10  cm  ,  the  device  will  behave  equivalently  to  many  smaller  capacitors  in 
parallel  each  with  a  slightly  different  V™.  This  will  cause  the  transition 
into  depletion  to  be  dispersed  along  the  voltage  axis  and  lead  to  an  incor¬ 
rectly  high  determination  of  interface  trap  density.  The  dotted  lines  in 
Figure  12  are  actual  quasistatic  and  1  MHz  data  which  are  compared  to  a 
calculated 


Figure  12.  Effect  of  surface  potential  nonuniformities  on  (.a)  quasistatic  C-V 
and  (b)  Calculated 


RPT50220 


16 


HUGHES 


SANTA  BARBARA  RESEARCH  CENTER 

1  SwflV  IM'V 

depletion)  and  in  some  of  these  cases  the  polarity  of  hysteresis  is  changed 
for  one  of  the  steps.  The  strong  flatband  shift  is  caused  by  charge  injection 
into  the  oxide  and  much  of  this  charge  is  mobile  as  witnessed  by  the  increased 
hysteresis.  A  tentative  explanation  for  the  steps  into  depletion  is  that 
charge  is  now  available  to  fill  and  empty  previously  unoccupied  oxide  states 
farther  from  the  interface.  This  could  also  account  for  the  occasional  shift 
in  hysteresis  polarity  by  contributing  mobile  ions  Cor  charges)  to  the  oxide. 

Monitoring  C-V  characteristics  versus  temperature  produced  results  that 
are  still  being  analyzed,  however  some  preliminary  conclusions  can  be  drawn. 
Figure  11  shows  a  1  MHz  C-V  curve  taken  at  77K  and  quasistatic  C-V  curves 
taken  at  four  temperatures  on  a  fairly  good  quality  MIS  structure 
(L)^t  *  1  x  10** cm-^  at  80K).  At  the  lower  temperatures  (<40K)  we  see  the 
widened  depletion  dip  and  the  inversion  spike  indicative  of  the  nonequilibrium 
measurement  of  an  interface  dominated  by  slow  insulator  traps.  This  occurs 


RAMP  VOLTAGE  (V) 

Figure  11.  1  MHz  and  quasistatic  at  7K.  Quasistatic  ramped  from  inversion  to 

accumraulation  versus  temperature  showing  discrete  states  at  onset 
of  inversion. 
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RAMP  VOLTAGE  (V) 

Figure  9.  Positive  bias  with  time.  V^g  follows  bias  indicative  of  charge 
injection  from  accumulation  layer. 


RAMP  VOLTAGE  (V) 

Figure  LU.  C-V  change  with  bias  during  cool  down.  Note  hysteresis  polarity 
change  in  4  to  6V  range  of  +5V  bias  and  "Steps"  in  depletion 
possibly  due  to  oxide  trap  occupation. 
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RAMP  VOLTAGE  (V) 

1  MHz  C-V  curves  on  Si  and  HCT.  Ramp  started  at  positive  bias. 

Si  hysteresis  opposed  to  bias,  HCT  hysteresis  follows  bias  indica¬ 
tive  of  mobile  ions  and  charge  injection,  respectively,  as  the 
major  contributors. 


Bias-stress  measurements  have  been  used  to  evaluate  the  stability  of  this 
charge  under  prolonged  bias  and  at  elevated  temperatures.  Standard  HgCdTe/ 
Photox”  Si02/Ti-Au  MIS  devices  were  tested  with  the  following  sequence: 
cooled  to  77K  at  0V  gate  bias;  biased  to  ±5V  and  ±10V  for  up  to  two  hours  with 
C-V  curves  taken  at  regular  intervals;  warmed  to  3U0K;  cooled  to  77K  with 
applied  biases  of  ±5V;  and  again  monitored  at  regular  intervals.  The  overall 
shift  seen  in  the  case  of  biasing  after  cool  down  was  in  the  direction  of 
applied  bias  and  occurred  mostly  in  the  first  five  to  ten  minutes,  as  the  data 
in  Figure  9  shows.  Very  small  shifts  opposed  to  the  applied  bias  were  some¬ 
times  seen  in  the  initial  two  or  three  minutes.  Thus  charge  injection  is 
still  predominant  under  prolonged  bias  stress  but  is  reasonably  stable  at  de¬ 
vice  operating  temperatures.  The  resulting  total  flatband  voltage  changed  only 
0.3V  in  two  hours.  Applying  bias  at  room  temperature  however,  produced  drama¬ 
tic  changes  in  the  C-V  curves  subsequently  taken  at  77K.  As  shown  in 
Figure  10,  flatband  voltages  were  shifted  strongly  (.4  to  5V ;  in  the  direction 
of  applied  bias,  the  magnitude  of  the  hysteresis  increased  a  factor  of  2  or  3, 
"steps"  appeared  (.apparently  indicative  of  an  incremented  transition  into  full 
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Figure  7.  Motion  of  charge  giving  rise  to  C-V  hysteresis,  (a)  positive 
charge  is  injected  into  oxide  requiring  stronger  negative  bias 
(b)  positive  charge  moves  closer  to  gate  requiring  less  negative 
bias. 

Upon  the  return  ramp  to  positive  bias,  these  charges  will  be  compensated  at  a 
less  negative  bias  than  was  previously  required  resulting  in  a  positive  flat- 
band  shift.  Charge  injection  on  the  other  hand,  will  result  in  a  negative 
flatband  shift  due  to  positive  charge  injection  into  a  previously  uncharged 
oxide  requiring  a  more  negative  bias  for  compensation.  This  logic  applies  to 
both  positive  and  negative  ions  and  biases.  The  result  is  that  C-V  hysteresis 
is  opposed  to  the  bias  for  mobile  ions  and  follows  the  bias  for  charge  injec¬ 
tion.  C-V  measurements,  with  the  voltage  ramp  initially  at  positive  bias,  are 
shown  in  Figure  8  for  both  silicon  and  HgCdTe.  Hysteresis  are  of  opposite 
polarities  indicating  mobile  ions  dominate  in  silicon  and  charge  injection 
dominates  in  HgCdTe.  The  amount  of  hysteresis  in  not  substantially  different. 
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(a)  ABRUPT  Si/SiO  INTERFACE  x  ALTERED  SURFACE 

L  REGION 

(b)  COMPLEX  HgCdTe  INTERFACE 


Figure  6.  Diagram  shows  comparison  of  abrupt  Si/SiO^  interface  with  HgCdTe/ 
SiC>2  interface  which  has  altered  (process-induced)  surface  region 
and  a  reduced  native  oxide  layer.  Origin  of  different  surface  and 
oxide  trapping  sites. 

Electrical  consequences  of  the  differences  between  the  two  systems  are 
apparent  in  trapping  characteristics  and  predominant  charge  distributions. 
Thermal  oxidation  of  the  silicon  creates  fixed  and  mobile  charge  in  the  oxide 
to  which  various  treatments,  such  as  HC1  during  oxide  growth,  have  been  ap¬ 
plied  as  remedies.  Mobile  charge  in  the  oxide  is  still  the  dominant  source  of 
instability.  Charge  injection,  however,  is  dominant  in  the  HgCdTe  system,  as 
can  be  seen  by  a  comparison  of  hysteresis  in  C-V  measurements.  (Hysteresis  is 
the  shift  in  the  C-V  curve  along  the  voltage  axis  after  biasing  positively  or 
negatively.  Figure  7  indicates  the  motion  of  charge  under  negative  bias  with 
(a)  charge  injection,  and  (b)  mobile  ions  as  the  source.  In  the  case  of  posi¬ 
tive  ions  in  Figure  7(b),  ions  be  moved  closer  to  the  negatively  biased  elec¬ 
trode. 
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spike  at  the  onset  of  inversion.  Both  effects  are  lessened  as  ramp  rate  is 
decreased,  indicating  that  equilibrium  has  not  been  reached.  It  is  of  course, 
possible  to  sit  at  each  voltage  until  equilibrium  is  obtained.  However,  use 
of  low  frequency  admittance  spectroscopy  will  give  a  more  complete  picture  of 
these  slow  states.  Whereas  C-V  measurements  use  a  fixed  frequency  and  scan  in 
voltage,  admittance  spectroscopy  selects  a  particular  voltage  and  scans  the 
frequency  range.  This  gives  the  technique  the  added  advantage  of  being  in 
equilibrium  and  enables  measurement  both  of  trap  (or  state)  distribution 
across  the  bandgap  as  well  as  unique  response  time  determination.  Admittance 

measurements  using  equipment  at  SRI  International,  capable  of  frequencies  as 

—4  —2 

low  as  10  Hz,  reveal  a  significant  trap  density  responding  in  the  10  Hz 

region. 

Interfaces  between  Photox™  and  HgCdTe  which  show  significant  concentra¬ 
tions  of  these  slow  states  have  been  measured  earlier  on  the  surface  study 
contract  and  reported  in  the  proposal  for  this  study.  These  are  interpreted 
to  be  states  in  the  insulator  which  communicate  with  the  HgCdTe  via  a  tunnel¬ 
ing  mechanism. 

Further  characterization  of  interface  traps,  oxide  states  and  their 
interactions  have  been  carried  out  in  bias-stress  measurements  and  C-V  versus 
temperature.  Bias-stress  measurements  will  indicate  sources  of  mobile  charge 
and  their  degree  of  mobility. 

It  is  instructive  to  compare  the  emerging  picture  of  the  Photox™/ HgCdTe 
interface  to  that  of  SiC^  on  Si.  On  a  physical  basis,  these  interfaces  differ 
dramatically.  A  schematic  of  the  two  systems  is  shown  in  Figure  b.  As  men¬ 
tioned  previously,  the  HgCdTe  surface  region  is  substantially  altered  stoi- 
chiometrically  from  the  bulk  by  polishes  and  etches  employed  in  surface  pre¬ 
paration.  This  surface  then,  on  exposure  to  air,  grows  a  native  oxide  »20A 
thick.  The  ensuing  passivation  of  the  surface  with  Photox™  Si02  exposes  the 
surface  initially  to  Silane  which  reduces  the  native  oxide  to  the  layer  of  Si, 

4 

Te  and  0  found  in  the  final  device.  Silicon,  on  the  other  hand,  is  a  sub¬ 
stantially  harder  material  with  only  one  constituent  and  hence  essentially 
maintains  bulk  crystallinity  up  to  the  surface.  When  thermal  oxidation  is 
performed,  a  clean  interface  between  the  silicon  and  the  SiC^  passivation 
results  with  no  intermediate  layer.  The  nature  of  the  SiC^  layers  are  also 
significantly  different  as  discussed  in  a  previous  report. 
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1  MHz  curves  showing  flatband  voltage  difference  between  n  and 
p  type  material. 
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Quasistatic  C-V  curve  with  ramp  rate  faster  than  required  for 
equilibrium.  The  widened  depletion  dip  gives  an  artificially  Low 
value  for  across  the  bandgap.  The  peak  at  the  onset  of 

inversion  will  not  occur  under  equilibrium  and  hence  makes  data 
reduction  here  meaningless. 
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Figure  3.  Measurement  of  dark  current  versus  gate  bias  for  gated  diode  (Gate 
overlaps  both  n  and  p  sides  showing  a  difference  in  the  voltage  at 
which  depletion  occurs  for  n  and  p  type  material.  After  J.P. 
Rosbeck  et  al. 

In  addition  to  the  fixed  charge  and  interface  trap  charge  discussed  a- 
bove,  the  existence  of  very  slow  states  has  been  detected  indirectly  through 
quasistatic  C-V  measurements  and  now,  seen  directly  with  admittance  spectro¬ 
scopy.  Measurement  of  D^t  is  typically  done  by  comparison  of  high  and  low 
frequency  C-V  measurements  and  use  of  Berglund  Integration  to  convert  gate 

bias  to  surface  potential  (¥  ).  This  requires  that  the  low  frequency  (or 

s 

quasistatic)  sufficiently  approximates  an  equilibrium  condition.  Typical  ramp 
rates  used  are  10-50  mV/sec  resulting  in  an  equivalent  frequency  on  the  order 
of  1  Hz.  Figure  5  shows  such  a  measurement  on  a  high  quality  MIS  device. 

Seen  here  are  an  artificially  wide  depletion  dip  and  a  capacitance 
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Figure  2.  Density  of  interface  traps  (Di{.)  derived  from  low  frequency  meas¬ 
urements  and  calculated  response  (ideal)  and  low  frequency  and 
high  frequency  measurements.  Superimposed  are  estimated  compo¬ 
nents  due  to  valence  and  conduction  band  tailing.  Shaded  areas 
show  charged  contributions. 


This  effect  has  been  seen  in  two  types  of  devices,  DLHJ  detectors  and  MIS  ca¬ 
pacitors.  Figure  3  shows  the  effect  on  dark  current  of  changes  in  gate  bias 
for  a  gated  photovoltaic  detector.  Dark  current  generation  from  interface 
states  will  increase  as  the  change  in  surface  potential  (via  changing  gate 
bias)  causes  depletion  on  the  p  or  n  side.  The  bias  at  which  this  occurs  for 
the  two  sides  of  the  junction  differs  by  nearly  1.3V  in  this  case,  as  seen  by 
the  two  peaks  in  dark  current.  Shown  in  Figure  4  are  capacitance-voltage 
characteristics  for  capacitors  fabricated  on  adjacent  wafer  sections  annealed 
to  comparable  doping  densities  of  n  and  p  type.  Again,  more  than  a  IV  differ¬ 
ence  is  seen  in  the  onset  of  depletion,  represented  here  as  the  flatband  volt¬ 
age  (VpB).  The  difference  in  Efs  also  causes  a  difference  in  the  work  func¬ 
tion  for  the  two  material  types.  However  this  effect  only  contributes  =>0.2V, 
leaving  the  major  contributor  to  be  the  different  interface  trap  charge. 
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Capture  cross  sections  were  measured  by  selecting  the  capacitor  bias, 
temperature,  and  rate  window  (  t)  which  gave  the  peak  position  corresponding  to 
the  trap  of  interest.  Gradually  shortening  the  fill  pulse  duration  ( 
enabled  measurement  of  the  trap  fill  time  as  the  slope  of  a  plot  of 
versus  signal  amplitude.*® 

RESULTS 

The  DLTS  spectra  showed  three  peaks  (Figure  1).  Two  were  very  prominent 
and  well  characterized  while  the  third  was  less  pronounced.  Obtaining  the 
activation  energies  of  the  trap  required  separating  the  capture  cross  section 
activation  energy  (Ea)  from  the  DLTS  value  (EqLts^  %LTS  =  Ef  +  Ea  where 

Et  is  the  actual  trap  activation  energy;  E a  may  be  determined  from  the  capture 
cross  section  values  and  fill  pulse  variation  data.  Trap  activation  energies 
were  determined  to  be  0.12  eV,  0.172  eV ,  and  0.079  eV.  These  are  electron 
traps  since  only  majority  carrier  traps  are  detected;  thus,  activation 
energies  quoted  are  measured  from  the  conduction  band  edge.  The  above 
energies  are  in  general  agreement  with  those  reported  for  traps  seen  in  bulk 
HgCdTe . 1 

The  activation  energies  for  the  two  prominent  peaks,  at  0.172  eV  and 
0.12  eV,  were  constant  at  lower  biases  and  then  decreased  as  bias  increased 
due  to  the  Poole-Frankel  effect.  The  0.12  eV  trap  activation  energy  dropped 
to  a  value  of  =0.05  eV  at  fields  of  10^V/cm.  The  activation  energies  quoted 
were  determined  from  low  bias  and  hence  are  independent  of  electric  field  so 
that  correction  for  the  Poole-Frankel  effects  is  not  required.  The  0.079  eV 
activation  energy,  is  also  reported  in  the  low  field  limit  and  is  assumed  to 
be  field  independent. 
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Selecting  the  experimental  parameters  listed  earlier  for  capture  cross 
section  measurement  resulted  in  a  capacitance  decay  curve  for  the  0.12  eV  peak 
which  has  two  slopes,  indicative  of  two  traps  emptying  at  different  rates. 
These  capture  cross  sections  are  3.5  x  10”*®  cm^  and  1.1  *  10“*®  cm^.  The 

same  method  applied  to  the  0.079  eV  peak  also  showed  two  traps.  One  capture 

cross  section  is  very  close  to  the  1.1  x  10“*®  cm^  obtained  previously  while 
the  second  showed  a  smaller  value  of  1.2  x  10  10  cm  .  As  can  be  seen  from  the 
DLTS  spectra  of  Figure  1,  there  is  some  overlap  of  the  peaks.  It  can  be 
inferred  therefore,  that  values  for  capture  cross  sections  obtained  from  the 
0.12  eV  peak  correspond  to  the  0.12  eV  and  0.172  eV  peaks  and  those  from  the 
0.172  eV  peak  correspond  to  the  0.172  eV  and  0.079  eV  peaks.  Due  to  the 
similarity  of  the  1.1  x  10  10  cm  value  obtained  in  both  measurements  this  was 
assigned  to  the  0.079  eV  peak.  This  leads  to  capture  cross  sections  for  the 
0.12  eV  and  0.079  eV  peaks  of  3.5  x  10“*®  cm^  and  1.2  x  10“*®  cm^, 
respectively.  (See  Table  1  for  summary.)  These  values  indicate  the  0.12  eV 

trap  to  be  neutral,  the  0.172  eV  trap  neutral  or  attractive  and  the  0.079  eV 


state  to  be  neutral  or  repulsive. 


The  capture  cross-sections  quoted  were 


measured  at  66K  and  102K  for  the  0.172  eV  and  0.12  eV  traps  respectively.  A 

n/KT 

temperature  dependence  of  o  =  '  is  assumed  in  this  region  but  has  not 

been  experimentally  verified. 

Choosing  a  bias  and  temperature  range  which  reveals  the  0.079  eV  peak  and 
plotting  the  capacitance  transient  logarithmically  showed  three  slopes  (Fig¬ 
ure  2),  revealing  three  distinct  decay  rates.  Plotting  a  family  of  curves 
with  temperature  as  the  parameter  showed  a  dependence  indicative  of  real  traps 
as  opposed  to  instrumental  effect.  The  third  segment,  however,  shows  an 
increased  slope  in  the  1,000  to  2,000  ys  region.  This  mandates  that  initia¬ 
tion  of  trap  decay  is  delayed  until  this  time.  Divalent  traps  may  explain 
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this  behavior,  with  one  trap  level's  decay  dependent  on  another's  emptying. 
Actual  details  of  the  trapping  mechanism  remain  to  be  determined. 

Depth  profiling  showed  the  0.172  eV  trap  to  have  a  fairly  constant 
concentration  from  the  maximum  depletion  width  up  to  the  surface.  The  0.12  eV 
trap  is  not  present  at  the  surface,  but  rises  in  concentration  to  a  level 
slightly  greater  than  that  for  the  0.172  eV  trap  at  0.5  un,  dropping  to  a 
nearly  equal  level  in  the  bulk  approximation  at  1.0  im,  as  shown  in  Figure 
3.  Due  to  the  very  low  relative  concentration  of  the  0.079  eV  trap,  as  well 
as  its  proximity  to  the  0.172  eV  trap,  an  accurate  profile  could  not  be 
determined. 

Since  DLTS  for  MIS  structres  simultaneously  profiles  two  trap  distri¬ 
butions;  one  in  energy  at  the  interface  and  the  other  spatially  from  the  inter¬ 
face.  Each  observed  peak  could,  in  general,  represent  a  bulk  trap  or  a 
surface  state.  To  differentiate  the  energy  from  the  spatial  distribution  we 
can  examine  the  energy  profile  for  surface  states  as  determined  by  C-V 
analysis.  This  U-shaped  energy  distribution  of  states  is  not  observed  in 
these  data,  either  as  a  set  of  discrete  states  or  as  two  state  distributions 
from  the  conduction  and  valence  bands.  A  single  discrete  surface  state,  on 
the  other  hand,  would  be  observed  only  when  the  band  bending  caused  the  Fermi 
level  to  cross  this  energy  level.  The  measured  activation  energies  of  0.12  eV 
and  0.172  eV  (again  measured  from  Ec) ,  and  the  biases  which  reveal  them,  are 
therefore  not  consistent  with  the  interpretation  of  surface  states;  these  are 
taken  to  be  trap  levels  within  the  bulk  material.  The  0.079  eV  peak  profiles 
were  not  complete  enough  to  be  able  to  def f erentiate  between  bulk  trap  and 
surface  state. 
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CONCLUSIONS 

For  the  first  time  we  have  been  able  to  examine  HgCdTe  electronic  proper¬ 
ties  near  its  interface  with  SiC^.  Three  distinct  trap  levels  have  been 
observed,  as  compared  to  only  two  seen  previously  in  the  bulk  material.  Bulk 
trap  levels  reported  *  were  at  l/2Eg  and  3/4Eg  with  capture  cross-sections  of 
3.5  x  10“^cm^  and  6.3  *  10”^cm^.  Within  the  limits  of  error  for  the  two 
sets  of  experiments  it  is  likely  that  these  are  the  same  traps  as  seen  in  this 
work  at  0.12  eV  and  0.172  eV.  The  non-uniform  distribution  with  depth  of  one 
trap  and  the  alteration  of  its  bulk  activation  energy  indicate  that  surface 
treatments  effect  bulk  trap  levels.  The  third  peak  detected  (0.079  eV)  may  be 
a  surface  state  to  which  DLTS  measurement  of  diodes  is  not  sensitive.  The 
information  obtainable  through  this  technique  represents  a  significant  in¬ 
crease  in  sensitivity  over  previously  reported  techniques  such  as  EER  and 
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Figure  Captions 


DLTS  Spectra 

Capacitance  Transient  with  Temperature  as  Parameter 
Relative  Trap  Concentrations  from  Depth  Profiling 
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Figure  2.  Capacitance  Transient  with  Temperature  as  Parameter 
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Section  II 
Channeled  Implants 

Channeled  implantation  of  HgCdTe  is  being  pursued  in  the  interest  of  pro¬ 
ducing  implanted  regions  where  the  doping  is  determined  by  the  implanted  ion 
concentration  and  not  by  the  resulting  lattice  damage.  Currently,  diodes  fab¬ 
ricated  using  ion  implantation  are  limited  to  n+  on  p  due  to  the  high  level  of 
n  type  damage  caused  by  the  nuclear  stopping  component  of  random  implants.  No 
satisfactory  method  of  annealing  out  this  damage  has  yet  been  developed. 

Atomic  profiles  of  channeled  implants  are  expected  to  have  a  random  peak 
at  the  surface  caused  by  ions  which  are  not  channeled  on  initial  impact  with 
the  surface,  followed  by  a  longer  "tail"  of  atoms  which  are  channeled  and 
hence  penetrate  more  deeply  into  the  material. 

Evaluation  of  these  profiles,  their  variation  with  mass  and  their  corres¬ 
ponding  electrical  activity  is  the  purpose  of  this  task.  The  techniques  used 
are  secondary  ion  mass  spectroscopy  (SIMS)  to  obtain  the  implanted  species 
atomic  profile,  and  capacitance-voltage  (C-V)  measurements  for  the  correspond¬ 
ing  electrical  activity.  Preliminary  evaluation  of  traps  created  in  the  dam¬ 
age  process  has  been  carried  out  using  DLTS. 

Samples  for  channeled  implant  are  grown  by  the  solid-state  recrystalliza¬ 
tion  method  (SSR),  which  produces  single  crystal  sections  of  HgCdTe.  These 
are  then  oriented  to  the  <110>  direction  and  sliced  into  wafers  parallel  to 
this  face.  Saw  damage  is  removed  by  a  Br/DMF  chemo-mechanical  polish.  An¬ 
nealing  at  240°C  in  Hg  vapor  for  four  weeks  results  in  a  doping  density  of 
1  x  lD^cm-^  n  type  to  a  depth  of  three  to  four  mils.  Several  ions  were 
looked  at  including  Na,  Cu  and  As.  Implant  energies  ranged  from  90  to  200  keV 
with  doses  of  1  *  lO^cm-^  to  5*  lO^cm”^.  Preparation  for  C-V  analysis  util¬ 
ized  the  standard  polish  and  etch  followed  by  an  additional  etch  to  remove  the 
desired  amount  of  material  for  depth  profiling.  MIS  devices  were  made  with 
Photox”  and  Ti/Au  gates. 

In  order  to  have  a  baseline  comparison,  random  implants  were  made  using 
unoriented  horizontal  zone  melt  material  at  similar  doses  and  energies.  These 
wafers  were  also  evaluated  using  SIMS  and  C-V. 
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The  STMS  profiles  of  channeled  implants  made  with  a  range  of  atomic 
masses  show  a  clear  trend  of  the  peak  position  moving  closer  to  the  surface 
with  increased  ion  mass.  This  is  shown  in  Figure  13.  For  the  same  dose  and 
energy,  the  amplitude  of  this  peak  also  increases  with  increasing  mass.  This 
is  the  expected  result  as  atomic  size  increases  relative  to  the  available 
channel  cross  section  more  ions  will  impact  the  lattice  and  be  implanted 
randomly.  Also  seen  is  the  corresponding  decrease  in  amplitude  of  the  chan¬ 
neled  "tail"  as  fewer  ions  are  available  for  channeling. 


DEPTH  (pm) 

Figure  13.  Atomic  density  profile  from  SIMs  showing 
change  in  distribution  with  atomic  size. 

A  comparison  of  random  and  channeled  SIMS  profiles  is  shown  in 
Figure  14(a)  for  Boron.  Here  the  peak  positions  are  at  the  same  depth,  with 
the  amplitude  somewhat  higher  for  the  random  case,  in  agreement  with  the 
assumption  that  this  portion  of  the  channeled  profile  represents  a  random 
component.  This  may  be  enhanced  in  our  case  by  the  polish  which  the  wafers 
receive  prior  to  implant  which  results  in  a  damaged  surface,  effectively 
plugging  the  lattice  channels. 
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Figure  14.  Comparison  of  electrical  activity  and  atomic  concentration 
profiles  (a)  for  Boron  (b)  for  Hydrogen. 


The  component  of  electrical  activity  measured  is  n  type  doping  concentra¬ 
tion  (Ny)  obtained  from  1  MHz  C-V  measurements.  At  this  frequency  interface 
traps  do  not  contribute  to  the  capacitance  in  accumulation  or  strong  inversion 
and  we  may  use 


N  4kT  e  2  C 

_ o _  _  ox  ox _ -2 

ln(N  /n.)  +  1/2  ln[21n(N  /n. )-l ]  q2e  X  2  Cu_(rain) 

Di  Di  soHF 


where  CQX  is  the  oxide  capacitance  measured  in  accumulation,  (min)  is  the 

minimum  capacitance  measured  in  strong  inversion,  xQ  is  the  oxide  thickness, 
eQ  and  eg  are  the  dielectric  functions  for  the  oxide  and  semiconductor  respec¬ 
tively  and  n^  is  the  intrinsic  carrier  concentration.  This  equation  assumes 
uniform  doping,  hence  in  our  case  it  yields  an  average  value  across  the  deple¬ 
tion  region.  This  width,  combined  with  depth  uniformity  and  measurement  er¬ 
rors,  are  reflected  in  the  error  bars  for  Figure  14  for  the  C-V  data.  The  er¬ 
ror  in  the  SIMS  measurement  is  =±1U%  of  the  final  sputter  depth,  about  U.23  pm 
in  this  case. 
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Electrical  activity  in  the  case  of  moderate  masses  such  as  boron  and 
phosphorous  appears  to  follow  the  random  component  of  the  profile  with  an 
activity  level  about  an  order  of  magnitude  higher  than  the  atomic  concentra¬ 
tion.  This  is  shown  for  boron  in  Figure  14(a).  The  implication  then  is  that 
to  first  order  that  electronic  stopping  does  not  contribute  any  significant 
electrical  activity.  Profiling  the  electrical  activity  of  random  implants 
appears  to  follow  a  similar  trend  with  a  slightly  shallower  distribution.  Due 
to  the  difference  in  doping  due  to  residual  impurity  levels  resulting  from  the 
two  growth  methods  (SSR  and  ZH)  for  these  samples,  higher  doses  were  used  for 
the  SSR  material  to  insure  electrical  activity  above  the  background  level. 


Figure  15.  Donor  concentration  due  to  implant  damage  scales 
linearly  with  0.23  eV  trap  concentration. 


it  is  possible  to  reliably  scale  atomic  concentration  for  changes  in  dose 
and  energy  using  LSS  theory.  The  mechanism  causing  damage  and  electrical 
activity,  however,  is  not  well  understood  rendering  extrapolation  impossible 
in  this  case.  Hence,  a  direct  comparison  of  doping  profiles  between  channeled 
and  random  implants  at  the  same  dose  and  energy  is  not  straightforward,  though 
some  data  can  be  scaled  for  comparison. 


% 


% 


% 


% 


RPT50220 


34 


di rect. 


All  of  the  prominent  structure  in  Figure  1  shows  some  dispersion  with  photon 
energy  in  the  normal  emission  geometry  and,  as  indicated  above,  can  reasonably  be 
accounted  for  in  the  bulk  band  structure  of  CdTe.  Thus,  the  measurements  are  not 
sensitive  to  a  state  of  distinctly  surface  origin.  The  lowest  binding  energy 
peak  in  He  I  (21.2  eV)  spectra  of  (110)  CdTe  has  been  previously  assigned  to 
surface  state  emission  [13,14].  The  observation  of  dispersion  of  this  peak  (peak 
6  in  Figure  1)  in  the  data  indicates  the  feature  has  a  substantial  bulk 
contribution.  At  the  lowest  photon  energies,  structure  6  has  shifted  far  from 

the  VBM;  a  surface  state  within  0.5  eV  of  the  VBM  may  account  for  the  nonzero 

emission  seen  in  the  curve,  but  the  possibility  of  optical  transitions  from 
near  the  band  edge  to  lower  lying  conduction  bands  must  also  be  examined. 

B.  HgCdTe 

The  result  of  assuminq  a  free  electron  final  state  band  with  G  =  il  (110) 

a 

and  the  same  value  of  '/0  as  was  used  for  CdTe  to  analyze  data  from  Hgo.ggWn. 3iTe 
is  shown  in  Figure  5  together  with  the  initial  states  derived  for  CdTe  in  this 
way  and  the  CdTe  bands  of  Ref.  22.  While  the  value  of  V0  appropriate  to  Hg-j_xCdxTe 
may  be  different  from  3.7  eV  in  the  final  analysis,  use  of  the  same  inner  potential 
facilitates  a  preliminary  examination  of  the  composition  dependence.  As  indicated 
in  Figure  2,  the  initial  state  structure  of  the  two  materials  is  highly  similar, 
with  features  in  the  alloy  occuring  at  increasingly  higher  binding  energy  as  the 
measurement  probes  further  from  the  VBM.  This  is  a  consequence  of  the  increasing 
s-electron  content  lower  in  the  bands  and  the  greater  binding  energy  of  the  Hg  6s 

level  compared  to  Cd.  When  the  data  from  the  x  =  0.31  sample  is  compared  to  that 

of  the  x  =  0.23  alloy  the  lowest  observed  state  (near  -5.4  eV)  and  bands  within 
3.0  e'/  of  the  VBM  are  found  to  be  fairly  insensitive  to  the  added  Hg  content 
while  the  structure  correspondi nq  to  the  one  labelled  3  in  Fig.  4  continues  to 
move  to  higher  binding  energy.  In  many  cases,  the  data  for  the  alloy  appears 
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bands,  other  values  of  G  may  contribute  [19].  Using  the  above  relation  of  k  and 

Ef  together  with  the  fact  that  the  intial  state  has  the  same  value  of  k  and  an 

energy  lower  than  Ef  by  hv ,  the  initial  state  bands  may  be  derived  from  the 

measured  final  state  energies  for  the  free  electron  parabola  with  G  =  £l[L  (HO). 

a 

These  results  are  given  in  Figure  4  together  with  the  calculated  valence  bands 
of  Chelikowsky  and  Cohen  [22].  The  value  of  V0  used  was  3.7  eV.  A  similar 
value,  4  eV,  has  been  obtained  from  an  analysis  of  the  primary  beam  energy  de¬ 
pendence  of  LEED  spot  intensities  [23]. 

The  correspondence  between  the  experimental  initial  states  and  the  calcu¬ 
lation  is  generally  good,  suggesting  the  theoretical  band  structure  correctly 
predicts  most  of  the  observed  features  in  the  ARPES  spectra.  Refinement  of  the 
analysis  may  further  improve  the  comparison.  The  shift  in  the  experimental 
points  in  structure  5  in  Figure  4  from  near  the  lower  of  the  two  split  bands  toward 
the  upper  one,  for  example,  results  from  a  change  in  the  relative  intensity  of 
emission  from  the  two  bands:  while  the  peak  position  was  selected  in  Figure  1, 
there  are  clearly  two  contributions  to  the  spectra  near  structure  5  between 
h  v=  16  and  18  eV  (Figure  1).  Proper  separation  of  these  components  would  give 
additional  detail  to  the  comparison  in  Figure  4.  The  most  important  correction 
to  the  analysis  to  produce  Figure  4  is  the  use  of  more  realistic  final  state 
bands.  The  structure  labeled  2,  the  portion  of  structure  1  from  k  =  0.5  to  0.7, 
and  structure  3  and  5  for  k  >  0.7  are  likely  the  result  of  transitions  to  final 
states  which  may  be  modeled  by  free  electron  parabolas  with  G  vectors  other  than 

ill  (110).  The  increase  in  intensity  in  the  oeax  labeled  3  in  Figure  1  as  the 
a 

final  state  changes  to  the  ii  (HO)  parabola  between  hv  =  17.5  and  21  eV 

a 

corresponding  to  0.6  _>  k  >_  0.35  in  Figure  4)  supports  this  explanation.  The 

di soersi onl ess  portion  of  structure  1  between  n.15  and  0.5  also  involves  final 

states  other  than  the  G  =  4 tt  (HO)  papabola,  but  these  transitions  may  not  be 

a 
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IV.  Discussion 


The  uncertainty  in  lo  [8]  renders  incomplete  the  determination  of  the 
quantum  numbers  for  electron  states  in  a  solid;  k ^ |  is  determined  by  the 
experimental  geometry,  but  k±  must  be  obtained  by  analysis.  If  calculated 
valence  bands  and  conduction  bands  of  sufficiently  high  energy  are  available, 
comparison  of  the  measured  initial  and  final  state  energy  with  theory  along  the 
appropriate  line  in  reciprocal  space  can  give  k_i_.  When  the  initial  state 
dispersion  only  is  known,  the  measured  peak  positions  can  be  used  to  map  the 
conduction  bands  [17,18].  To  derive  the  valence  hand  E  versus  k  relation  when 
the  final  state  structure  is  also  not  known,  on  the  other  hand,  approximate 
conduction  bands  such  as  free  electron  parabolas  may  prove  successful  [19], 
Electronic  structure  calculations  for  CdTe  in  the  (110)  direction  have  been 
performed;  [20-22]  however,  final  states  of  sufficient  energy  are  not  available. 
The  present  ARPES  data  can  be  compared  with  calculation  if  a  free  electron  final 
state  is  assumed.  By  such  a  procedure,  the  accuracy  of  the  initial  state  calcu¬ 
lations  can  be  studied,  and  deviations  from  a  free  electron  parabola  for  the 
high  lying  conduction  bands  examined.  Preliminary  results  for  the  alloy  samples 
can  be  obtained  from  similar  analysis. 

A.  CdTe 


In  the  case  where  the  final  states  lie  on  a  free  electron  parabola,  Ef(k)  is 


given  by 


Ef  .-|_(k*G)2-V0 


where  G  is  a  reciprocal  lattice  vector  and  V0 ,  the  energy  difference  between  the 

bottom  of  the  free  electron  parabola  and  the  valence  band  maximum,  is  the  inner 

potential  [8],  The  strongest  contribution  to  the  photocurrent  should  come  from 

the  final  state  band  with  G  =  0  in  the  extended  zone  scheme  or  the  portion  of  the 

parabola  translated  to  the  first  Brillouin  zone  by  G  =  (110)  [19]  for  the 

a 

range  of  photon  energies  used.  Because  of  the  complexity  of  the  true  final  state 
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The  similarity  observed,  however,  suggests  that  complete  comprehension  for  CdTe 
or  an  alloy  of  a  given  composition  will  greatly  aid  analysis  for  alloys  with 
different  x  value. 

In  addition  to  location  of  the  electron  states  in  a  solid  by  means  of  analysis 
of  the  position  and  movement  of  peaks  in  the  spectra,  information  as  to  the 
symmetry  of  the  states  can  be  obtained  from  the  polarization  dependence  of  the 
peak  intensity.  Figure  3,  spectra  of  CdTe  and  HgQ# 77(^0. j^Te  for  the  three 
polarizations  used  and  hv  =  16  eV,  illustrates  this  sensitivity.  (The  data 
has  been  normalized  to  the  incident  flux  and  corrected  for  differing  backgrounds 
in  the  region  of  interest.)  For  s-polarization  along  the  [0011  direction,  one 
observes  considerable  enhancement  of  the  peak  near  -2.7  eV  in  both  cases  but 
particularly  in  CdTe,  while  the  structure  closest  to  the  VBM  is  not  strongly 
excited  for  this  polarization.  When  oriented  in  the  [001]  direction,  the 

polarization  vector  is  even  upon  reflection  in  the  (T 1 0)  mirror  plane  of  the 

% 

crystal.  Since  in  normal  emission  the  final  state  wavefunction  is  similarly 
even,  a  nonzero  matrix  element  for  optical  excitation  requires  the  initial  state  to 
also  have  even  reflection  symmetry  [12].  Thus  the  valence  band  state  contributing 
near  -2.7  eV  must  be  even  upon  reflection  in  the  (110)  plane  while  the  states 
contributing  nearest  the  VBM  at  this  photon  energy  are  of  odd  character.  While 
the  trends  observed  with  changing  polarization  reflect  the  overall  similarity 
of  the  electronic  structure  of  the  alloy  for  all  compositions  studied,  the  detailed 
behavior  depends  on  the  initial  states  involved  for  a  given  photon  energy  and  the 
wavelength  dependent  optical  properties  of  each  material,  which  determine  the 
electric  field  direction  inside  the  crystal  and  thus  the  excitation  [8],  The 
symmetry  of  the  states  contributing  to  the  ARPES  measurement  as  determined  by  the 
polarization  dependence  of  the  spectra  can  be  used  as  an  important  test  of  the 
comparison  of  the  experimental  electronic  structure  to  theory. 


previously  reported  normal  emission  data  from  (110)  CdTe  [13,14]  recorded  with 
the  He  I  resonance  line  (21.2  eV)  indicated  reasonable  agreement,  particularly 
given  differences  in  the  polarization  of  the  exciting  radiation.  Dispersion  of 
the  observed  structure  with  the  exception  of  the  highest  binding  energy  peak  near 
-4.6  eV  was  similarly  reported  in  ARPES  studies  of  fixed  hv  and  variable  emission 
angle  [14]  as  well  as  in  angle  integrated  experiments  with  variable  hv  [5]. 

The  dispersion  with  photon  energy  of  features  in  the  electron  energy  distri¬ 
bution  curves  from  CdTe  seen  in  Figure  1  is  characteri Stic  of  all  samples  studied. 
Although  the  general  features  in  the  spectra  from  HgTe  are  the  same  as  found  for 
CdTe,  the  dispersion  is  somewhat  less  pronounced  most  likely  as  a  result  of  the 
rougher  surface  produced  by  a  poorer  cleave  [15].  (Time  did  not  permit  additional 
cleaves  to  be  made  in  an  effort  to  improve  the  surface  quality.)  Indeed,  a  corres¬ 
pondence  of  the  structure  observed  for  CdTe  and  the  alloys  in  terms  of  relative 
location  and  polarization  and  photon  energy  dependence  can  be  found  in  each  case. 

A  comparison  of  the  emission  along  the  sample  normal  for  p-pol arization  of  hv  = 

18  eV  for  the  four  samples  studied  is  given  in  Figure  2.  For  the  crystals  containing 
Hg,  the  highest  binding  energy  feature  is  found  near  -5.5  eV  independent  of 
composition,  (If  the  virtual  crystal  description  of  the  electronic  structure 
held  for  these  states,  a  monotonic  change  in  energy  would  be  expected  [15].)  The 
intense  peak  in  the  curve  for  CdTe  near  -3.5  eV  is  observed  at  smoothly  increasing 
binding  energy  in  the  alloys  and  HgTe.  Similarly  the  two  peaks  or  peak  and 
shoulder  near  the  VBM  are  replicated  for  each  composition,  while  structure  between 
-2  and  -3.5  eV  apparently  increases  in  strength  as  the  Hg  content  increases. 

Because  the  energy  and  intensity  of  features  in  the  data  depend  on  the  final 
state  energy  and  dispersion  as  well  as  that  of  the  initial  states,  analysis  of 
the  final  state  bands  for  each  compound  must  be  performed  in  order  to  fully 
understand  the  composition  dependence  of  the  curves  at  a  given  photon  energy. 
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and  constant  wavelength  resolution  AA  =  8.7  A.  Samples  of  CdTe  (II-VI,  Inc.) 
and  solid  state  recrystall  ized  Hg-|_xCdxTe  (x  =  0.00,  0.23,  and  0.31)  were 
introduced  into  the  vacuum  chamber  following  bakeout  through  an  air  lock.  Clean 
(110)  surfaces  of  the  5  x  5x  ~  1 5  mm^  bars  were  exposed  by  cleaving  across  the 
square  faces  in  vacuum  (p  ~  2  x  10*^  Torr).  The  cleaved  surfaces  of  the  alloy 
were  found  to  be  n-type  by  determination  of  the  Fermi  level  position.  Departure 
from  the  bulk  stoichiometry  was  not  resolved  in  UPS  spectra ;  the  ratio  of  Hg  5d 
to  Cd  4d  emission  in  the  data  for  the  alloys  scaled  with  bulk  composition  in  the 
expected  way.  Prior  to  ARPES  measurements,  the  crystal  axes  were  oriented 
using  the  Low  Energy  Electron  Diffraction  (LEED)  pattern. 

III.  Results 

Electron  energy  distribution  curves  for  normal  emission  from  (110)  CdTe 
under  p-polarized  irradiation  (26°  from  normal  along  [110])  are  presented  in 
Figure  1  as  a  function  of  initial  state  energy  relative  to  the  VBM  for  0.5  eV 
increments  of  the  photon  energy.  Prominant.  peaks  in  the  data  are  marked  and  are 
observed  to  shift  in  binding  energy  and  modulate  strongly  in  intensity  as  the 
photon  energy  is  varied.  Such  peaks  occur  when  valence  band  and  conduction  band 
states  of  the  same  value  of  k  in  the  reduced  Brillouin  zone  scheme  are  separated 
in  energy  by  the  photon  energy.  As  hv  is  varied,  the  values  of  k  and  initial 
state  energy  for  which  these  conditions  are  met  change  in  a  way  characteristic 
of  the  dispersion  of  the  initial  and  final  state  bands  involved.  An  exception 
to  this  type  of  behavior  is  seen  in  the  peak  near  -4.6  eV  (structure  1)  at  the 
higher  photon  energies.  This  apparently  stationary  peak  arises  from  the  mainly 
cation  s-electron  derived  states  which  have  little  dispersion  deep  in  the  valence 
bands  and  near  the  zone  boundary.  The  characteri Stic  motion  of  the  peaks  marked 
in  Figure  1  will  be  used  to  compare  the  measured  electronic  structure  to  calculation 
in  Section  IV.  Comparison  of  spectra  for  photon  energies  of  21.0  or  21.5  eV  to 


% 


the  symmetry  line  through  the  points  TKX  in  the  Brillouin  zone  of  the  face 
centered  cubic  lattice  [11].  A  second  geometry  was  used  to  restrict  emission 
from  the  mainly  s-electron  derived  state  3.5  to  6  eV  below  the  VBM  (where  the 
nonperiodic  potential  effects  are  strongest)  to  the  hexagonal  face  of  the  8rillouin 
zone  along  the  WLW  line  [11].  This  data,  which  bears  on  the  breakdown  at  the 
virtual  crystal  approximtation  at  the  zone  boundary,  will  appear  in  a  later 
publ ication. 

Measurements  in  these  geometries  were  carried  out  in  an  ADES  400  (VG 
Scientific)  angle-resolving  system.  The  u-metal  chamber  is  pumped  by  a  400  1 /s 
triode  ion  pump  and  titanium  sublimation  pumping.  Electron  energy  analysis  is 
accomplished  by  a  moveable  hemispherical  electrostatic  analyzer  with  channeltron 
electron  multiplier  detector  and  pulse  counting  electronics.  The  analyzer  accep¬ 
tance  cone,  2°  degree  in  half  angle,  can  be  freely  oriented  relative  to  the 
sample  by  rotating  the  analyzer  from  5°  above  to  90°  below  the  horizontal  plane 
and  by  rotating  the  analyzer  through  360°  about  the  vertical  chamber  axis. 

Operating  with  the  pass  energy  fixed  provided  a  constant  energy  resolution  of 
0.2  eV. 

Linearly  polarized  radiation  with  1 3  _<  hv  £  27  eV  from  the  8°  beam 
line  (1-1)  at  the  Stanford  Synchrotron  Radiation  Laboratory  served  to  excite  the 
photoelectrons.  By  varying  the  angle  of  incidence  of  the  photon  beam  and  rotating 
the  sample  about  the  surface  normal,  the  direction  of  polarization  could  be 
controlled.  The  orientation  of  the  polarization  influences  the  spectrum  through 
symmetry  dependent  optical  selection  rules  [12]  governing  the  transition  from 
initial  to  final  state.  Measurements  were  made  in  the  normal  emission  geometry 
with  increments  in  the  photon  energy  of  0.5  or  1.0  eV  for  the  vector  potential 
oriented  70°  from  the  sample  normal  (stronger  s-pol arization)  along  both  the 
[001]  and  [110]  directions  and  26°  from  the  normal  (stronger  p-pol arization) 
toward  [110].  The  choice  of  monochromator  slit  settings  gave  a  2  x  2  mm2  spot 
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bond  [6].  Based  on  evidence  of  the  binding  energy  shifts  of  the  Hg  5d  and  Cd  4d 
core  levels  between  the  alloy  and  the  binary  compounds,  the  disparity  in  the 
bond  strength  of  the  two  metal  atoms  increases  upon  alloy  formation,  with 
important  implications  in  areas  such  as  heterojunction  formation  and  dislocation 
creation  and  motion  [7]. 

To  probe  the  dependence  of  the  electronic  structure  on  composition  and 
explore  the  effects  of  the  disorder  in  the  cation  potential  on  the  valence 
bands  in  the  alloy,  angle-resolved  photoemission  spectroscopy  (ARPES)  [8] 
is  a  particularly  well-suited  technique.  The  anisotropy  of  the  allowed  electron 
states  of  a  solid  in  reciprocal  space  is  reflected  in  the  spatial  variation  of 
the  photoexcited  electron  emission.  By  measurement  of  the  energy  distribution 
of  the  electrons  photoemitted  into  a  small  solid  angle  in  a  particular  direction 
relative  to  axes  in  the  crystal,  information  on  both  the  energy  and  crystal  momen¬ 
tum,  k,  of  the  valence  and  conduction  band  states  in  the  solid  can  be  obtained. 
ARPES  has  been  used  in  this  way  to  map  the  valence  and  high  lying  conduction 
band  states  in  a  range  of  materials  [9]  including  metal  alloys  such  as  CuNi  [10], 
which,  like  HgCdTe,  exhibits  strong  potential  disorder.  In  the  present  study, 
ARPES  has  been  used  to  examine  the  electronic  structure  mainly  along  the  (110) 
symmetry  line  of  CdTe  and  HgTe  as  well  as  Hg^_xCdxTe  alloys  with  x  =  0.23  and  x  = 
0.31.  The  spectra  obtained  are  presented  in  section  III  and  a  comparison  of 
the  valence  electronic  structure  of  the  materials  is  given  in  section  IV. 

II.  Experimental 

The  electron  states  to  which  one  is  sensitive  in  an  angle-resolved  experiment 
are  a  function  of  the  photon  beam  energy  and  polarization  and  the  position  of 
the  electron  detector.  The  spectra  reported  here  were  recorded  with  the  electron 
analyzer  along  the  normal  to  the  (110)  cleavage  face  of  the  samples  studied.  In 
this  configuration ,  kjj  is  required  to  he  zero  and  all  initial  states  lie  along 
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I.  Introduction 


Accurate  knowledge  of  the  band  structure  of  a  material  provides  the  means 
to  estimate  properties  ranging  from  the  ease  of  dislocation  formation  to  the 
conduction  band  effective  mass.  The  HgTe-CdTe  superlattice  system  is  a  case  in 
point:  electronic  [1]  and  optical  [2]  properties  of  the  material  have  been 
predicted  theoretically  prior  to  observation  which  may  make  the  superlattice  of 
technological  interest.  The  accuracy  of  the  quantitites  derived  from  the 
electronic  structure  depends  on  how  well  the  band  structure  is  known,  motivating 
study  of  the  electron  states  in  a  given  material. 

The  relation  between  the  electronic  structure  and  materials  properties 
is  particularly  evident  in  Hgi_xCdxTe.  The  strong  relativistic  interactions 
between  the  nucleus  and  the  valence  s-electrons  in  atomic  Hg  coupled  with  the 
influence  of  the  just  completed  d-electron  shell  serve  to  increase  the  binding 
energy  of  the  Hg  6s  electrons  by  1.4  eV  [3]  compared  to  the  lighter  Cd,  where 

these  effects  are  smaller.  As  a  result  of  the  shift  in  the  valence  electron 

binding  energy,  the  conduction  band  minimum  in  HgCdTe,  which  is  cation  s-electron 
derived,  moves  closer  to  the  valence  band  maximum  (VBM)  as  Hg  is  substituted 
for  Cd  in  the  lattice  [4].  At  the  same  time,  the  VBM  is  comprised  of  Te  p-derived 

states  and  is  fairly  insensitive  to  the  composition  change.  Thus  the  ability  to 

tune  the  alloy  bandgap  by  varying  the  composition  follows  from  the  valence  elec¬ 
tronic  properties  of  the  constituents  [5,6,7], 

A  second  immediate  consequence  of  the  alloy  electronic  structure  relates 
to  the  bonding  of  Cd  and  Hg  in  the  lattice.  The  disorder  in  the  potential  in 
the  cation  sublattice  introduced  by  the  differing  strength  of  the  interaction 
of  the  s-electrons  with  the  cation  nuclei  results  in  distinct  contributions 
to  the  alloy  valence  electronic  structure  from  Hg  and  Cd  [5].  Consequently, 
the  Cd  bond  to  Te  in  the  alloy  differs  from  that  of  Hg.  In  the  respective 
binaries,  the  Hg-Te  bond  is  considerably  more  metallic  and  weaker  than  the  Cd-Te 
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ABSTRACT 


Angle-resolved  photoemission  spectroscopy ,  a  technique  which  provides  direct 
information  on  the  quantum  numbers  k  and  E  which  characterize  the  electron  states 
in  a  periodic  solid,  has  been  used  to  study  the  electronic  structure  of  the  alloy 
Hgi_xCdxTe  for  x  =1.0,  0.31,  0.23  and  0.0.  Oata  were  collected  mainly  for 
emission  along  the  normal  to  the  cleaved  (110)  face  of  these  crystals  as  a 
function  of  photon  energy  with  various  directions  of  polarization.  The  structure 
observed  in  the  spectra  of  the  alloy  are  similar  in  dispersion  and  polarization 
dependence  to  that  obtained  for  CdTe.  Preliminary  analysis  of  the  observed 
features  using  a  free  electron  final  state  reveals  the  bulk  origin  of  the 
measured  transitions  by  comparison  to  calculated  valence  bands  and  elucidates 
aspects  of  the  composition  dependence  of  the  valence  band  states  involved  in 
the  optical  excitation. 
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The  unusual  bond  length  relaxation  exhibited  in  HgCdTe  is  discussed 
in  section  III.  The  weak  Hg-Te  bond  leads  to  mechanical  instability  of 
the  HgCdTe  lattice;  the  trends  predicted  for  the  bond  length  relaxation 
in  this  alloy  should  exacerbate  this  instability.  In  order  to  test  our 
overall  theoretical  approach,  the  theoretical  model  presented  in  this 
section  is  shown  to  predict  the  trend  observed  experimentally  by  Mikkelson 
and  Boyce  [3]  for  the  pseudobinary  alloy  InAs-GaAs.  This  work  is  included 
as  an  appendix  to  section  III. 

The  effects  of  this  unusual  bonding  behavior  on  the  thermodynamics 
and  kinetics  governing  growth  at  interfaces  are  seen  experimentally  in 
the  last  two  sections.  The  heterojunction  between  LPE  HgCdTe  is 
studied  in  the  fourth  section,  where  it  is  seen  that  HgCdTe  exhibits  a 
more  metallic-like  recovery  behavior  than  CdTe  when  bombarded  with 
ions.  The  observations  also  support  the  weakness  of  the  Hg-Te  bond 
and  the  consequent  ease  of  formation  of  a  vacancy  attributed  to  Hg. 

The  kinetics  of  the  activated  oxygen  uptake  on  the  surface  of  the 
HgCdTe  are  investigated  in  section  V.  The  uptake  of  oxygen  for  the 
first  two  monolayers  appears  to  scale  with  the  number  of  Hg-Te  bonds 
in  the  alloy.  Hence,  the  breaking  of  this  bond  and  the  subsequent 
removal  of  Hg  from  the  lattice  appears  to  play  an  important  role 
in  the  beginning  stages  of  oxidation. 
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I.  Introduction 

The  underlying  theme  brought  forth  in  our  fundamental  studies  of 
the  HgCdTe  surface,  interface,  and  bulk  regions  is  the  unusual  nature 
of  the  bonding  in  the  alloy.  The  differences  in  the  cation  bonding 
[1]  in  the  alloy  are  manifest  in  the  electronic  structure  and  mechan¬ 
ical  behavior  of  the  material.  Theoretical  calculations  [2]  reveal 
that  not  only  does  the  Hg-Te  bond  exhibit  more  metal  1 ic -1 i ke  behavior 
than  is  typically  observed  in  tetrahedrally  bonded  covalent  semicon¬ 
ductor  systems,  but  that  Hg-Te  and  Cd-Te  bond  length  relaxation  in 
the  alloy  exhibits  unusual  trends  as  well.  The  trend  is  for  the  already 
shorter  Hg-Te  bond  to  shorten  as  Hg  is  added  to  CdTe  and  the  longer 
Cd-Te  bond  to  lengthen  as  Cd  is  added  to  HgTe  .  This  suggests  that 
an  unusual  micro-mechanical  behavior  exists  in  the  alloy  that  may 
markedly  affect  the  bulk  thermodynamic  properties  of  the  material. 

This  should  therefore  produce  significant  effects  in  growth  processes. 

For  example,  those  that  govern  materials  growth  at  LPE  interfaces 
as  well  as  in  bulk  regions. 

Studying  the  electronic  structure  of  HgCdTe  yields  valuable  insights 
into  bonding  nature  and  electronic  behavior.  Once  the  bar,J  tructure 
is  understood,  many  properties  of  the  alloy  may  be  calculated.  Section 
II  discusses  the  experimental  determination  and  our  current  knowledge 
of  the  electronic  structure  of  HgCdTe. 
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In  the  case  of  hydrogen,  Figure  14(b),  the  electrical  activity  is  dis¬ 
tinctly  n  type  and  well  above  the  background  level  due  to  residual  impurities, 
however  there  is  no  apparent  correlation  between  atomic  and  electrical  pro¬ 
files  of  channeled  samples.  SIMS  of  random  Hydrogen  implants  is  not  yet  a- 
vailable. 

DLTS  measurements  on  samples  implanted  with  Cl  at  a  range  of  doping  den¬ 
sities  (due  to  implant  damage)  show  a  new  trap  with  an  activation  energy  of 

0.20  eV  (E  =  0.27  eV  at  200K  where  the  DLTS  peak  is  found)  measured  from  the 

°  -20  2 
conduction  band  edge.  The  capture  cross  section  of  9.6  x  10  cm  is  indica¬ 
tive  of  a  weak  repulsive  trap.  As  shown  in  Fig  16  the  concentration  of  this 
trap  decreases  almost  linearly  as  doping  density  decreases,  which  is  to  be 
expected  if  it  is  due  to  damage  created  by  the  implant.  The  concentration  of 
this  trap  is  an  order  of  magnitude  lower  than  the  doping  density  which  gives 
it  a  nearly  one-to-one  correspondence  with  implant  species  concentration. 

_  i  o  o 

Another  trap  is  also  seen  with  E.  =  0.053  eV  and  o  =  2  x  10  cm  which 
becomes  more  dominant  as  doping  density  approaches  the  background  level.  This 
trap  is  most  likely  the  same  as  the  0.079  eV  level  seen  in  unimplanted 
material.  Neither  of  the  prominent  traps  as  seen  in  unimplanted  material  at 
0.12  eV  and  0.172  eV  have  been  detected  in  these  implanted  samples. 
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to  show  less  dispersion  than  observed  for  CdTe  (see  Figure  5);  this  nay  be  an 
effect  of  the  non-periodic  alloy  potential. 

As  is  evident  in  Figure  5,  all  peaks  in  the  data  from  Hg q . 61  Cdf). 31  Te  (and 
the  other  samples  as  well)  can  be  accounted  for  by  the  bulk  rather  than  surface 
electronic  structure  in  analogy  to  the  assignments  of  these  features  to  the 
valence  bands  of  CdTe.  Of  the  published  band  structures  for  HgCdTe,  those 
which  give  the  valence  bands  in  the  (110)  direction  [20,24,25]  predict  band 
widths  which  are  generally  less  than  observed.  Recent  calculations  performed  in 
the  coherent  potential  approximation  [26,27]  are  the  most  suitable  for  comparison 
to  the  experimental  data;  however,  the  valence  bands  in  the  (110)  direction  are 
not  given. 

V.  Summary 

The  utility  of  ARPES  lies  in  the  sensitivity  of  this  technique  to  the 
quantum  numbers  k  and  E  which  characterize  the  electron  states  of  a  periodic 
solid.  ARPES  data  collected  normal  to  the  (’ •  0)  face  of  CdTe  and  HgCdTe  alloys 
have  been  examined  by  assuming  a  free  electron  final  state.  The  true  final  states 
in  these  materials  are  clearly  more  complex,  as  illustrated  by  the  divergence  of 
experiment  and  a  calculated  band  structure  for  some  features  in  the  data.  Still, 
much  of  the  observed  initial  state  dispersion  corresponds  to  prediction,  tentatively 
confirming  the  accuracy  of  the  calculation  as  well  as  the  nearly  free  electron 
nature  of  some  final  states.  Refinement  of  the  comparison  to  theory  by 
i ncorporati on  of  additional  final  states  and  consideration  of  the  photon  beam 
polarization  dependence  will  allow  greater  accuracy  in  the  determination  of  the 
correctness  of  calculation  and,  ultimately,  of  predicted  materials  parameters 
derived  from  the  electronic  structure. 
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Figure  Captions  . 


I  Figure  1  ARPES  spectra  recorded  along  the  normal  to  the  (110)  face  of 

CdTe  with  1  3  <_  hv  <_  27  eV  in  0.5  eV  steps.  The  radiation  was 
incident  at  64°  to  the  normal  (stronger  p-pol arization)  with 
1  the  electric  field  direction  along  [T 1 0] . 


Figure  2 


I 


Comparison  of  normal  emission  spectra  (p-polarization)  taken  at 
hv  =  18  eV  for  the  four  Hg-|_xCdxTe  compounds  studied.  The 
curves  show  highly  similar  structure  which  tends  toward  higher 
binding  energy  with  increasing  Hg  content. 


Figure  3 


Figure  4 


The  symmetry  of  the  initial  states  is  reflected  in  the  polarization 
dependence  of  the  structure  in  the  normal  emission  spectra  at  a 
given  photon  energy  (16  eV).  The  enhancement  of  the  peak  near 
-2.7  eV  and  suppression  of  emission  near  the  VBM  with  hv 
polarized  parallel  to  [001]  is  observed  in  all  samples  (CdTe  and 
H00.77CdO.23Te  are  shown.) 

Initial  state  dispersion  relations  for  the  normal  emission  data  from 
CdTe  derived  assuming  a  free  electron  final  state  (points)  are  com¬ 
pared  to  the  calculated  bands  of  Ref.  22  (solid)  in  the  (110) 
direction.  See  text  for  discussion  of  comparison. 


Figure  5  Initial  state  dispersion  relations  derived  from  the  experimental 

features  for  CdTe  (open  points)  and  Hgo.69CdQ.31Te  (filled  points) 
and,  as  a  guide,  the  calculated  bands  for  CdTe  (Ref.  22).  The 
addition  of  Hg  increases  the  binding  energy  of  most  features;  and 
the  alloy  potential  disorder  may  reduce  the  measured  dispersion 
of  some  peaks. 
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Section  1 1 1  - 1 1 1 


BOND-LENGTH  RELAXATION  IN  HgT_xCdxTe  ALLOY 

The  alloy  system  Hg]_xCdxTe  (MCT)  exhibits  unusual  bond-length  relaxation 
behavior.  In  most  of  the  pseudo-binary  alloys  in  which  the  bond-length  differ¬ 
ence  of  the  two  components  is  fairly  large,  such  as  Ga-|_xInxAs  or  ZnSei_xTex, 
the  behavior  of  the  bond-length  relaxation  can  be  predicted  by  a  simple  radial 
force  model  (Appendix  1).  In  that  model  the  bond-length  relaxation  is  deter¬ 
mined  by  minimizing  the  distortion  energy  which  results  from  the  deviation 
of  the  bonds  from  their  "natural  bond-lengths",  where  the  "natural  bond-lengths" 
are  defined  as  the  bond-lengths  of  the  pure  components.  However,  in  the  case 
of  MCT  in  which  the  bond-lengths  of  the  two  components  are  nearly  the  same 
(within  0.3%),  the  effect  considered  in  the  previous  model  is  small;  the 
dominating  factor  comes  from  the  coupling  of  the  bonding  state  and  the  neighbor¬ 
ing  anti-bonding  states,  which  is  termed  "metallization"  in  Harrison's  bonding 
theory  [1], 

Harrison  [2]  showed  that  in  calculating  the  bond  energy  of  HgTe  and  CdTe, 
the  metallization  term  which  is  neglected  in  the  Bond  Orbital  Approximation  is 
no  longer  small  and  cannot  be  neglected.  In  the  MCT  alloy  the  bonds  have 
different  local  environments,  and  the  metallization  terms  change  depending  on 
the  neighboring  atoms.  We  will  show  that  the  changes  of  the  metallization 
terms  are  responsible  for  the  behavior  of  the  bond-length  relaxation  in  MCT. 

In  Harrison's  theory  [1]  the  bond  tension  T  is  defined  as  the  first  deriv¬ 
ative  with  respect  to  internuclear  distance  of  all  terms  contributing  to  the 
bond  energy  excluding  the  overlap  interaction.  The  bond-length  is  determined 
by  the  balance  of  the  tension  and  the  repulsion  force. 

Now  let  us  consider  an  alloy  Hg-j.xCdxTe  with  x  near  1.  This  is  a  dilute 
solution  of  Hg-Te  in  a  CdTe  host  (Figure  la).  We  will  compare  this  with  the 
Hg-Te  bond  in  pure  HgTe,  (Figure  lb)  and  investigate  the  effect  of  the  variation 
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of  the  metallization  term  with  environment.  The  metallization  terms  contributing 
to  the  bond  tension  of  the  Hg-Te  come  from  the  couplings  between  the  Hg-Te 
bonding  state  and  its  neighboring  anti-bonding  states;  and  the  couplings  between 
the  neighboring  bonding  states  and  the  Hq-Te  anti-bonding  state.  In  the  alloy 
(Figure  1),  the  couplings  through  the  Hg  cation  are  identical  to  the  case  in 
the  pure  HgTe,  and  the  couplings  through  the  Te  anions  are  different  from  the 
pure  HgTe  case.  In  pure  HgTe  the  tension  terms  on  the  Hg-Te  bond  which  arise 
from  the  couplings  through  the  Te  anion  are:  3Tman-  [3],  which  is  the  coupling 
between  the  Hg-Te  bonding  state  and  the  three  neighboring  Hg-Te  antibonding 
states;  and  3Tmaj/i  which  is  the  coupling  between  the  neighboring  Hg-Te 
bonding  states  and  the  Hg-Te  anti-bonding  state.  Whereas  in  the  MCT  the  tension 
terms  on  the  Hg-Te  bond  which  arise  through  the  Te  anion  are:  3TrT)ai2'  >  The 
couplings  between  the  Hg-Te  bonding  state  and  the  neighboring  Cd-Te  anti-bonding 
states;  and  2Tma]'2  which  is  the  coupling  between  the  neighboring  Cd-Te 
bonding  states  and  the  Hg-Te  antibonding  state.  Therefore,  the  change  in  the 
tension  on  the  Hg-Te  bond  in  MCT  with  respect  to  the  pure  HgTe  is 
A  Tm  =  3  (Tmai2"  -  Tmall'  +  Tmal'2  "  Tmal'l) 

These  tension  terms  can  be  Obtained  by  taking  the  first  derivative  of  the 
coupling  energies  with  respect  to  the  Hg-Te  bond  length.  These  are  given  by: 
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where  apn  and  acn  are  defined  as  the  polarity  and  covalency  for  bond  n; 
ean  and  ebn  are  energy  levels  of  anti-bonding  state  and  bonding  states  of 
bond  n;  V2  and  V3  are  covalent  energy  and  polar  energy  of  Hg-Te;  and  V-|_  is  the 
metallic  energy  of  Te. 

When  we  examine  these  formulae,  we  find  that  the  behavior  of  i  Tm  is  pre¬ 
dominantly  determined  by  the  ap's.  The  polarity  of  Hg-Te  is  smaller  than 
that  of  Cd-Te,  i.e.  a  pi  <  a  p2*  One  then  finds  Tman'  <  Tmai  2'  and 

Tmal'l  <  ^mal'2  hence  A  Tm>  0  for  Hg-Te  in  the  alloy.  Similarly  we 

will  find  that  A  Tm  <  0  for  Cd-Te  in  the  alloy.  A  positive  A  Tm  for  Hg-Te  in 

the  alloy  predicts  the  bond-length  will  decrease  from  that  in  the  pure  Hg-Te. 

The  opposite  trend  is  predicted  for  Cd-Te  in  the  alloy. 

To  calculate  a  numerical  value  of  A  Tm,  precise  atomic  term  values  are 
needed.  Once  the  value  of  A  Tm  is  obtained,  the  final  bond-length  can  be 
obtained  in  one  of  the  following  ways: 

(1)  Use  the  experimental  spring  constant  K  of  the  bond.  In  this  case 

A  d  =  -  A  Tm/k. 

(2)  Plug  the  A  Tm  back  into  the  balance  equation  Vo(d)  =  Td/4  [1],  where 
Vo(d)  is  the  overlap  repulsion  energy  and  proportional  to  d"4. 

(3)  Use  numerical  methods  to  find  the  bond  energy  v.s.  interatomic  distance. 
The  bond-length  and  bond  energy  in  the  alloy  can  thus  be  obtained  simul 
taneousl y . 

The  first  method  is  only  semiempi rical .  The  second  method  does  not  require  the 


experimental  spring  constant,  and  it  results  in  a  sel f-consi stent  theory.  How- 


is  required,  and  besides  being  computationally  more  tedious,  only  yields  correct 
trends  rather  than  accurate  values.  In  the  third  method,  expressions  for  each 
term  contributing  to  the  bond  energy  and  use  of  a  computer  are  required.  On 
the  other  hand,  this  technique  yields  other  useful  bonding  properties.  Sher  et 
al  [4]  have  used  the  third  method  and  obtained  the  result  that  for  the  case  of 
dilute  Hg-Te  in  the  Cd-Te  host,  the  bond-length  of  Hg-Te  decreases  by  0.02  A  ; 
for  dilute  CdTe  in  the  HgTe  host,  the  bond-length  of  Cd-Te  increases  by  0.02  A  . 

The  calculated  numerical  value  of  the  bond-length  variation  suffers  from 
difficulties  in  obtaining  precise  atomic  term  values.  However,  by  understanding 
the  physical  terms  which  are  responsible  for  the  bond-length  relaxation,  we  can 
predict  the  trend  in  MCT.  The  important  result  is  that  the  shorter  Hg-Te  bonds 
become  shorter  and  the  already  longer  Cd-Te  bonds  are  lengthened.  This  trend 
is  opposite  to  that  observed  in  III-V  systems  such  as  InGaAs,  where  the  effects 
of  metallization  can  be  ignored.  Future  work  will  concentrate  on  the  implica¬ 
tions  of  the  bond-length  relaxation  on  the  lattice  stability  and  the  bonding  of 
a  heterojunction  interface. 

References 

1.  W.  A.  Harrison,  Phys.  Rev.  B,  27,  6  (1983)  3592. 

2.  W.  A.  Harrison,  J.  Vac.  Sci.  Technol.  A  _1_,  (1983)  1672. 

3.  In  Tman  'the  'm'  stands  for  this  tension  come  from  the  metallization  the 
'a'  stands  for  the  coupling  through  the  anion. 

4.  A.  Sher,  A-B  Chen,  W.  E.  Spicer,  C-K  Shih,  Proceedings  of  the  1984  U.S. 
Workshop  on  the  Physics  and  Chemistry  of  Mercury  Cadmium  Telluride  to  appear 
in  0.  Vac.  Sci.  Technol.,  Feb.,  1985. 


is  shown  in  Figure  2a,  where  a  high  density  of’ extrinsic  stacking  faults 
(ESF)  is  evident.  Examination  of  the  HgCdTe  side  of  the  heterojunction 
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reveals  a  significant  difference  in  damage  behavior  when  compared  to  that 
of  the  CdTe.  At  this  higher  magnification,  the  defects  in  the  HgCdTe 
were  much  harder  to  locate  than  in  the  CdTe.  Very  few  defects  are  exhibited 
on  this  side  of  the  heterojunction,  where  a  representative  area  of  the 
HgCdTe  reveals  no  defects  as  shown  in  Figure  2b.  The  defects  studied  in 
the  HgCdTe  generally  had  little  or  no  dissociation  of  the  1/6<112>  Shockley 
partials.  Figure  2c  shows  an  almost  perfect  60°  dislocation  in  the  HgCdTe, 
where  the  arrow  points  down  the  extra  {111}  half-plane.  To  the  best 
of  our  knowledge,  perfect  60°  dislocations  have  not  been  observed  in  CdTe. 

ESF  are  also  observed  in  the  HgCdTe,  a  typical  fault  is  shown  in  Figure 
2d,  but  the  separation  between  the  two  Shockley  partial  dislocations  is 
always  much  smaller  in  this  material  than  is  observed  in  CdTe. 

IV.  Discussion 

These  results  obtained  from  the  TEM  examination  of  ion  beam  milled 
HgCdTe  and  CdTE  are  unexpected  and,  at  first  glance,  surprising.  Because  of 
its  mechanical  instability  compared  to  CdTe,  one  might  predict  that  the 
HgCdTe  would  exhibit  higher  defect  densities  than  CdTe  after  being  subject  to 
identical  milling  conditions.  In  addition,  the  CdTe  stacking  fault  energy 
(calculated  to  be  =■  lOmJ/m?  [14])  is  predicated  to  be  greater  than  that  of 
HgCdTe  on  the  basis  of  the  usual  effect  of  alloying  on  stacking  fault  energy 
in  metals  [15].  That  there  is  a  difference  in  the  cation-Te  bonding  in  the 
CdTe  and  HgCdTe  as  is  indicated  by  experimental  [7]  and  theoretical  consider¬ 
ations  [16,17]  is  supported  by  these  observations. 

The  extremely  high  defect  density  in  the  ion  milled  CdTe  is  not  sur¬ 
prising  since  this  magnitude  of  damage  is  typically  observed  by  Yamashita 


general  distribution  of  native  and  induced  defects  in  the  two  materials 
was  examined  at  low  magnification  using  conventional  2-beam  imaging. 

Since  the  HgCdTe  and  CdTe  are  ion  milled  under  identical  conditions 
during  the  final  thinning  step,  a  direct  comparison  of  the  density, 
distribution  and  nature  of  the  ion  bombardment  induced  defects  can  be  made 
between  HgCdTe  and  CdTe.  In  addition,  the  HgCdTe/CdTe  heterojunction  can 
be  examined.  Figure  1  shows  typical  defect  structures  in  the  HgCdTe, 

CdTe,  and  interfacial  region,  respectively.  These  micrographs  reveal 
a  quite  surprising  result:  under  identical  ion  beam  milling  conditions, 
the  CdTe  contains  a  much  higher  density  of  defects  (2-3  orders  of  magnitude) 
than  the  adjacent  HgCdTe.  Figure  la  reveals  that  the  defects  in  the  HgCdTe  side 
of  the  heterojunction  consist  of  single-line  dislocations  (located  at  A) 
and  occasional  dislocation  loops  (located  at  B) .  The  defect  density  increases 
over  several  microns  to  the  heterojunction  where  it  abruptly  increases  in  the 
interfacial  region,  as  seen  in  Figure  lc.  The  long  line  dislocations  are  seen 
to  pile  up  at  the  interface  on  the  HgCdTe  side  of  the  metallurgical  junction, 
and  the  damage  induced  defect  loops  increase  in  density  from  the  HgCdTe  towards 
the  junction.  In  the  interfacial  region,  these  loops  greatly  increase  in 
density  and  become  smaller.  This  high  density  stays  approximately  constant 
into  the  CdTe  substrate,  as  shown  in  Figure  lb. 

Invaluable  insight  into  the  nature  of  the  ion  beam  induced  defects  can 
be  obtained  from  high  resolution  images  taken  from  these  regions.  Lattice 
images  of  the  heterojunction  region  taken  normal  to  the  {110}  plane 
and  looking  down  the  {111}  planes  are  shown  in  Figure  2.  Inspection 
of  the  CdTe  side  of  the  heterojunction  reveals  a  very  high  density  of 
dissociated  60°  dislocations  and  interstitial  loops,  consistent  with  the 
above  observations  on  their  density.  A  representative  area  of  the  CdTe 


nique  should  be  feasible  for  HgCdTe  if  low  ion  beam  energies  are  employed 
and  the  sample  is  cooled  to  liquid  nitrogen  temperatures  [11].  With  such 
TEM  samples,  investigations  into  native  defect  distribution  and  elemental 
composition  profiles  in  epitaxial  HgCdTe  grown  on  CdTe  and  alternate 
substrates  may  be  performed.  It  is  hoped  that  additional  bonding  information 
can  be  obtained  in  these  studies.  This  paper  will  report  on  a  preliminary 
TEM  investigation  supplementing  our  ongoing  study  of  the  bonding  nature 
of  HgCdTe. 

II.  EXPERIMENTAL 

Cross-sectional  TEM  (XTEM)  specimens  were  prepared  from  a  sample  of 
HgCdTe  (x  =  0.3)  grown  by  LPE  on  a  {111}  CdTe  substrate.  The  LPE 
material  was  grown  from  a  Te  melt  to  a  layer  thickness  of  14.3  pm. 

Most  of  the  material  used  in  this  study  was  supplied  by  Santa  Barbara 
Research  Center  (SBRC);  however,  studies  in  Rockwell  samples  give  similar 
results.  Employing  a  XTEM  specimen  technique  similar  to  that  of  Yamashita 
on  samples  of  thin  films  grown  on  CdTe  [12],  XTEM  specimens  suitable  for 
viewing  down  the  plane  of  the  heterojunction  interface  were  prepared. 

The  samples  were  oriented  such  that  a  <110>  direction  lay  normal  to  the 
plane  of  sample,  so  that  little  tilting  was  needed  to  align  the  electron 
beam  with  the  <1 1 0>  zone  axis  for  high  resolution  lattice  imaging  [13]. 

Final  thinning  was  achieved  by  ion  milling  at  room  temperature  to  perforation 
using  Ar+  at  4  keV  with  a  35  pA-mm“2  beam  current  incident  at  15°  for 
several  hours.  Examination  of  the  samples  was  made  using  a  Philips 
EM400ST  having  a  high  resolution  objective  lens. 

III.  Results 

Using  the  Energy  dispersive  analysis  of  X-rays  (EDX)  facility  the 
interfacial  region  between  the  HgCdTe  and  CdTe  was  located,  and  the 
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defect  properties  exhibited  in  this  LPE  heterojunction  suggest  that  an 
unusual  interface  kinetics  and/or  chemistry  exist  in  this  system.  In 

light  of  the  results  mentioned  above,  a  question  may  be  raised  as  to 

J 

the  importance  of  lattice  match  at  the  HgCdTe/CdTe  interface.  Instead, 
questions  of  chemical  driving  forces  across  the  interface  may  be  critical. 

For  example,  the  unusual  bonding  of  HgCdTe  provides  a  possible  driving 
force  for  moving  Cd  from  CdTe  into  HgCdTe  [7,8],  Many  scientific  clues 
into  the  mechanisms  of  the  bonding  in  the  alloy  may  exist  in  this  hetero¬ 
junction.  By  studying  this  heterojunction  in  depth,  it  is  hoped  that  a 
more  fundamental  understanding  into  the  nature  of  the  bonding  in  this 
alloy  may  be  obtained. 

Transmission  electron  microscopy  (TEM)  is  a  useful  tool  for  direct 
defect  structural  visualization  and  X-ray  compositional  analysis.  Prepa¬ 
ration  of  TEM  specimens  requires  thinning  techniques  such  as  ion  beam 
milling.  Ion  milling  is  particularly  useful  on  the  preparation  of  cross- 
sectional  samples,  which  are  useful  in  studying  defect  structures  at 
interfaces  or  defect  distributions  associated  with  ion  implanted  material. 
While  cross-sectional  TEM  has  been  successfully  employed  to  great  effect 
in  the  elemental  [9]  and  III-V  [10]  semiconductors,  difficulties  have 
been  encountered  in  implementing  this  technique  to  the  HgCdTe  system  due 
to  the  damage  introduced  into  the  crystals  during  ion  milling.  This 
induced  damage  tends  to  mask  the  original  defect  structure  and  hence  impede 
native  defect  structural  analysis.  Insight  into  the  nature  of  the  bonding 
in  the  alloy  might  be  gained,  however,  if  the  ion  milling  defects  are 
studied. 

Recently,  progress  in  refining  the  ion  milling  technique  for  application 
to  HgCdTe  has  been  made,  and  preliminary  reports  indicate  that  this  tech- 
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I.  INTRODUCTION 

In  order  to  exploit  the  attractive  electronic  capabilities  of  the 
Hg-|_xCdxTe  alloy  system,  detailed  defect  structural  characterization  of 
this  material  must  be  undertaken.  As  outlined  by  Cheung  [1]  in  these 
proceedings,  numerous  studies  have  shown  that  the  electronic  behavior  of 
the  alloy  is  linked  to  the  mechanical  defect  structure  introduced  during 
the  growth  and  processing  of  the  material.  Potential  device  performance 
is  often  impeded  by  the  presence  of  this  defect  structure.  Of  extreme 
interest  are  those  defects  associated  with  interfaces,  and  in  particular, 
the  defect  structure  that  arises  at  the  liquid  phase  epitaxial  (LPE) 

HgCdTe  grown  on  CdTe  heterojunction. 

Some  of  the  compositional  and  structural  properties  that  characterize 
this  heterojunction  are  unique  among  LPE  heterojunction  systems.  The 
transition  region  width  between  the  HgCdTe  and  the  CdTe  substrate  has  been 
shown  to  be  primarily  a  function  of  the  quality  of  the  CdTe  surface  preparation; 
but  even  the  best  quality  LPE  HgCdTe/CdTe  heterojunctions  are  characterized 
by  transition  widths  on  the  order  of  one  micron  [2,3].  What  are  the  fundamental 
mechanisms  leading  to  such  a  great  transitional  width?  For  comparison,  typical 
LPE  heterojunction  widths,  such  as  the  AlGaAs/OaAs  system,  are  on  the  order  of 
100  \  [4],  The  lattice  mismatch  in  that  system  is  roughly  one-third  of  that 
on  the  HgCdTe/CdTe  system,  but  the  LPE  transition  width  is  approximately  100 
times  less  that  that  in  the  LPE  HgCdTe/CdTe.  In  addition,  analysis  of 
the  HgCdTe/CdTe  LPE  heterojunction  shows  an  integrated  misfit  dislocation 
intensity  in  the  so-called  misfit  zone  that  is  much  greater  than  can  be 
explained  by  lattice  mismatch  between  the  HgCdTe  and  CdTe  [2,5].  These 
misfit  dislocations  are  not  completely  confined  to  the  {111}  growth  plane 
in  the  HgCdTe  side  of  the  heterojunction  [6].  The  unique  structural 
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TEM  INVESTIGATION  OF  THE  OIFFERENCES  IN  ION  MILLING  INDUCED  DAMAGE 
OF  Hg-j _xCdxTe  AND  CdTe  HETEROJUNCTIONS 
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Abstract 

The  lattice  defects  induced  by  ion  milling  at  the  liquid  phase  epitaxy 
(LPE)  grown  HgCdTe  on  CdTe  heterojunction  have  been  investigated  using 
transmission  electron  microscopy  (TEM).  TEM  analysis  reveals  that  subjecting 
both  materials  to  identical  milling  conditions  results  in  much  higher  defect 
densities  (2-3  orders  of  magnitude)  in  the  CdTe  than  in  the  HgCdTe  (x  =  0.30). 
In  addition,  high  resolution  lattice  images  of  both  regions  reveal  that  the 
stacking  fault  widths  in  the  HgCdTe  are  much  narrower  than  in  the  CdTe,  which 
implies  that  the  stacking  fault  energy  is  much  higher  in  the  HgCdTe  as  com¬ 
pared  to  the  CdTe.  These  observed  phenomena  are  shown  to  support  the 
weakness  and  metallicity  of  the  Hg-Te  bond  as  proposed  by  theoretical  con¬ 
siderations. 
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Figure  Captions 

Figure  1  (a)  An  Indium  atom  substitutes  for  Gallium  in  GaAs;  the  tetrahedral 
symmetry  is  preserved  at  the  In  atom. 


(b)  The  displacement  of  As  atoms,  u,  results  in  the  Ga-As  bonds 
being  decreased  by  u  cos  9.  The  In-As  bond  lengths  are  varied 
dGaAs  +  u  -  dInAs 
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By  minimizing  the  energy  with  respect  to  u  we  obtain 
u  =  3/4(dInAs  -  dGaAs), 

or  the  deviation  of  the  In-As  bond-length 
AdInAs  =  1/A  (dGaAS  "  dInAs)* 

Similarly  the  bond  length  deviation  of  dilute  Ga-As  in  an  InAs  host  is 
AdGaAs  =  1/4  (dInAs  “  dGaAd)« 

This  predicts  exactly  Mikkelsen  and  Boyce's  result  [1].  Note  that 
this  displacement  leaves  all  metallic  atoms  in  their  original  sites;  only 
the  As  atoms  move.  To  a  good  approximation  this  is  true  in  the  concen¬ 
trated  alloy  where  the  metal  atoms  form  a  face  center  cubic  lattice  of 
the  weighteci  average  lattice  parameter  and  the  As  atoms  move  to  fit  local 
bond  lengths. 

One  can  apply  this  principle  to  the  rocksalt  structure  systems;  e.g. 
Ki-xRbxBr,  and  predict  that  the  bond-length  deviation  of  the  solute  in  the 
dilute  solution  is  half  of  the  difference  between  the  two  components. 

A  dKBr  =  l/2(dRbBr  -  d«Br) 

A  dRbBr  =  l/2(dKBr  -  dRbBr) 

However,  because  of  the  lattice  geometry,  the  approximation  that  all  the 
cations  stay  in  their  original  sites  might  not  be  as  good.  In  rocksalt 
structure,  the  restoring  force  on  the  displaced  first  neighbor  results 
from  the  atom  right  behind  the  first  heighbor;  one  would  expect  that  the 
relaxation  extends  more  than  one  interatomic  distance  and  the  value  of 
the  bond-length  deviation  would  be  somewhat  smaller  than  what  is  predicted 
by  this  model.  This  agrees  with  the  experimental  result  [2]. 
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In  calculating  pseudo-binary  alloy  properties,  the  nearest  neighbor 
distance  is  a  very  important  parameter.  For  many  years  people  have  used  the 
Virtual  Crystal  Approximation  (VCA)  which  assumes  that  all  atoms  occupy  the 
average  lattice  positions.  In  1982  Mikkelson  and  Boyce  used  Extended  X-ray 
Absorption  Fine  Structure  (EXAFS)  to  study  the  atomic  scale  structure  of  the 
Gai.xMxAs  system  [1].  They  found  that  the  Ga-As  and  In-As  bond  lengths  did 
not  follow  the  VCA;  instead,  they  only  slightly  deviated  from  their  natural 
bond  lengths  even  in  the  dilute  solutions,  and  the  magnitude  of  the  bond- 
length  deviation  at  the  dilute  limit  was  about  a  quarter  of  the  difference  of 
the  natural  bond  lengths  of  the  two  components  (GaAs  and  InAs).  These  results 
can  be  explained  by  a  simple  radial  force  constant  model  in  which  the  weaker 
angular  forces  are  neglected. 

Consider  a  compound  GaAs  in  which  one  of  the  Ga  atoms  is  replaced  by  an 
In  atom.  This  configurationTepresents  the  dilute  limit  of  In-As  in  a  GaAs  host 
(see  Figure  1).  In  this  configuration,  four  In-As  bonds  are  identical  and  the 
tetrahedral  symmetry  is  preserved  at  the  In  atom.  We  regard  any  difference 
between  Ga  and  In  as  a  first  order  difference.  Let  us  call  dQa^s  the  natural 
bond  length  of  the  host  and  K  the  stretching  force  constant  for  the  host. 

Adding  In  will  change  K  by  an  amount  we  regard  as  first  order,  but  that  alone 
will  not  give  the  distortion.  However,  because  the  natural  bond  length  of 
InAs,  din/\s,  is  larger  than  d qa/\s  by  a  first  order  amount  the  As  atoms  around 
the  In  impurity  move  outward  by  a  first  order  distance;  we  hold  the  second 
neighbors  fixed  and  to  first  order  can  take  K  to  be  the  value  for  the  host. 

Let  us  call  u  the  outward  relaxation  displacement  of  the  As:  Then  Ad jn As ,  the 
amount  that  the  final  In-As  bond  length  deviates  from  dinAS,  is  dQa/\s  -  djn/\s 
+  u.  The  bond  lengths  of  the  neighboring  Ga-As  are  decreased  by  -u  cos  0 
equal  to  -u/3  to  first  order.  The  energy  associated  with  this  distortion  is 


BOND-LENGTH  RELAXATION  IN  PSEUDO-BINARY  ALLOYSt 

C.  K.  Shi h 
W.  E.  Spicer* 

W.  A.  Harrison 

Stanford  University,  Stanford,  CA  94305 
Arden  Sher 

SRI  International,  Menlo  Park,  CA  94025 

ABSTRACT 

The  bond-length  relaxation  in  pseudo-binary  alloys  can  be  predicted  by 
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Figure  Captions 


Figure  la,b  Dark  circles  represent  Te  atons,  circles  with  dots  in  them 
represent  Hg  atoms,  and  open  circles  represent  Cd  atoms. 
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in  work  done  on  CdTe  [12],  but  the  relatively  low  resulting  damage  to  the 
HgCdTe  lattice  is  quite  unexpected.  In  fact,  this  low  defect  density  is 
suggestive  of  a  metallic-like  recovery  behavior  in  the  HgCdTe.  Cole  et 
al .  [18]  showed  that  the  room  temperature  hardness  of  HgCdTe  (x  «*  0.3) 
and  CdTe  are  similar,  suggesting  that  the  rate  of  glide  of  dislocations  is 
similar  in  the  two  materials.  Therefore  the  differences  in  the  observed 
defect  structures  on  either  side  of  the  junction  is  not  attributed  to  the 
differences  in  the  glide  behavior. 

The  difference  in  the  distribution  of  defects  in  this  ion  milled 
HgCdTe/CdTe  heterojunction  may  be  understood  by  considering  the  mechanism 
of  climb  in  these  two  materials.  The  lattice  images  from  both  of  these 
regions  reveal  that  there  is  a  smaller  dissociation  of  the  1 /6<1 1 2> 
Shockley  partial  dislocations  in  the  HgCdTe  relative  to  the  CdTe.  Since 
climb  occurs  by  a  diffusion-controlled  process  in  which  vacancies  or 
constituents  must  diffuse  to  the  dislocation  core,  the  climb  of  the 
narrower  stacking  faults  and  perfect  60°  dislocation  in  the  HgCdTe  should 
be  much  easier.  This  implies  that  the  rate  of  recovery  would  be  much 
higher  in  the  HgCdTe  than  in  the  CdTe,  as  is  observed. 

This  recovery  mechanism  might  also  ^e  aiued  by  the  large  density  of 
point  defects  present  in  the  ion  bombarded  HgCdTe  lattice  arising  from  the 
weak  Hg-Te  bond  [17]  and  the  consequent  ease  of  formation  of  the  cation 
vacancy  attributed  to  the  weakly  bound  Hg  [8,16].  Hence  this  "Hg  vacancy" 
mechanism  is  consistent  with  taese  theoretical  calculations.  Thus  it  is 
possible  that  both  materials  undergo  comparable  damage  during  ion  milling, 
but  that  the  presence  of  a  high  density  of  induced  "Hg  vacancies"  along 
witn  toe  narrower  stacking  faults  in  the  HgCdTe  enhance  the  recovery  rate 
in  this  material  over  that  of  CdTe.  One  must  also  ask  whether  the  ease 
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of  movement  of  dislocations  in  HgCdTe  may  contribute  to  the  extreme 
(>  1  urn)  heterojunction  width  in  the  HgCdTe  regions  of  the  HgCdTe/CdTe 
heterojunctions. 

V.  SUMMARY 

Preliminary  analysis  of  the  ion  bombarded  induced  defects  in  HgCdTe/CdTe 
heterojunctions  using  XTEM  yield  interesting  phenomena  which  manifest  the 
differences  in  the  cation-Te  bonding  occurring  in  HgCdTe  and  CdTe.  Under 
identical  ion  milling  conditions,  the  CdTe  contains  much  higher  resulting 
damage  (2-3  orders  of  magnitude  more  defects)  in  the  form  of  small  inter¬ 
stitial  or  vacancy  loops  compared  to  HgCdTe.  Lattice  images  reveal  that 
the  separation  of  the  Shockley  partials  is  in  general  much  less  in  the 
HgCdTe  than  in  the  CdTe.  These  observations  suggest  that  the  HgCdTe  under¬ 
goes  a  much  higher  rate  of  recovery  than  the  adjacent  CdTe,  where  the 
narrower  stacking  faults  would  tend  to  quicken  the  rate  of  climb  of  the 
defects  and  consequently  lead  to  an  increased  rate  of  annihilation.  This 
mechanism  would  be  enhanced  by  the  high  density  of  point  defects  in  the 
form  of  "Hg  vacancies"  that  would  be  present  in  the  ion  bombarded  HgCdTe 
due  to  the  weak  Hg-Te  bond.  Hence,  these  observations  support  the  theoretical 
work  which  predicts  a  weak  Hg-Te  bond  and  consequent  ease  of  formation  of 
the  "Hg  vacancy"  in  the  alloy. 
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Figure  Captions 

Figure  1(a)  Conventional  2-beam  image  of  the  ion-milled  HgCdTe,  showing 
line  defects  (A)  that  were  believed  to  be  present  in  the  LPE 
material,  and  small  dislocation  loops  (B)  that  were  intro¬ 
duced  during  ion  milling. 

(b)  TEM  micrograph  of  the  CdTe  substrate  side  adjacent  to  the 
heterojunction,  showing  large  densities  of  small  loops,  and 
long-line  native  defects. 

(c)  2-beam  image  of  the  interfacial  region  clearly  showing  the 
difference  in  ion  beam  milling  damage  in  both  sides  of  the 
metallurgical  junction. 

Figure  2(a)  TEM  high  resolution  lattice  image  taken  on  the  CdTe  side  of 
the  heterojunction.  This  representative  area  is  seen  to 
possess  a  high  density  of  extrinsic  stacking  faults  (ESF). 

(b)  Lattice  image  taken  from  a  representative  area  on  the  HgCdTe 
side  of  the  heterojunction  showing  no  defects. 

(c)  Lattice  image  showing  an  almost  perfect  60°  edge  dislocation 
in  the  HgCdTe  side  of  the  heterojunction.  The  arrow  points 
down  the  extra  half-plane. 

(d)  Lattice  image  slowing  ESF  in  the  HgCdTe.  The  widths  of  the 
ESF  in  the  HgCdTe  were  generally  much  less  than  in  the  CdTe. 
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COMPARATIVE  STUDY  OF  ACTIVATED  OXYGEN  UPTAKE  ON  HgTe,  Hg0.69Cd0  31Te, 

AND  CdTe 
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ABSTRACT 

Soft-X-ray  photoemission  spectroscopy  was  used  to  monitor  the  uptake 
of  oxygen  on  the  (1.10)  faces  of  HgTe,  Hgo.59Cdo.3i  Te,  and  CdTe  during  the 
formation  of  roughly  two  monolayers  of  oxide.  The  oxidation  was  activated 
by  a  hot  filament  in  1  i re-of-sight  of  the  sample  surface.  The  oxides  formed 
contain  Cd  and  Te;  whereas,  Hg  is  lost  from  the  surface  during  the  initial 
stages  of  oxide  formation.  The  rate  of  oxidation  for  the  Hg0.6gCdQ.3i  Te 
sample  was  found  to  be  69%  as  great  as  that  for  the  binary  HgTe  and  lower 
than  observed  for  CdTe.  The  scaling  of  oxidation  rate  with  the  number  of 
Hg-Te  bonds  between  HgTe  and  the  x  =0.31  alloy  together  with  the  observed 
release  of  Hg  upon  oxidation  suggests  the  importance  of  breaking  the 
Hg-Te  bond  as  a  step  in  the  oxidation  process.  Evidence  for  an  unstable 
character  of  the  oxide  on  the  Hg  containing  compounds  is  also  presented. 


a)  Fannie  and  John  Hertz  Foundation  Fellow 

b)  Present  address:  Ministry  of  Defense,  P.  0.  Box  2250,  Haifa,  Israel 

c)  Stanford  W.  Ascherman  Professor  of  Engineering 


88 


% 


I.  Introduction 

In  addition  to  its  relevance  to  passivation  processes,  study  of  the 
oxidation  of  Hg]_xCdxTe  can  afford  insights  into  the  manifestation  of  the 
alloy  electronic  structure  in  the  chemical  activity  of  the  surface.  Experi¬ 
mental  [1]  and  theoretical  [2]  studies  of  the  alloy  electronic  structure 
indicate  distinct  Hg  and  Cd  contributions  to  the  valence  states  4  to  6  eV 
below  the  valence  band  maximum,  implying  a  difference  in  the  bonding  of  Cd 
and  Hg  in  the  lattice.  The  bond  in  HgTe  is  considerably  more  metallic  and 
weaker  than  the  bond  in  CdTe  [3],  Upon  alloy  formation,  a  net  charge 
transfer  from  Cd  to  Hg  takes  place,  as  evidenced  by  shifts  in  the  binding 
energy  of  the  Hg  5d  and  Cd  4d  core  level  photoemission  peaks  [4]  and  chemical 
shifts  measured  as  a  function  of  composition  in  nuclear  magnetic  resonance 
data  [5].  Such  a  charge  rearrangement  results  in  a  higher  electron  density 
near  the  Hg  and  a  stronger  Cd-Te  bond  in  the  alloy  at  the  expense  of  the 
bond  of  Hg  to  Te  [4,6].  These  features  of  the  alloy  electronic  structure 
may  be  expected  to  influence  the  oxygen  chemistry  of  the  surface  through 
modification  of  the  kinetics  of  reactions  which  occur  or  by  determining  the 
nature  of  the  oxide  formed.  A  comparative  investigation  of  the  activated 
uptake  of  oxygen  on  atomically  clean  (110)  surfaces  fo  CdTe,  HgTe,  and 
**90.69^0.31  Te  was  undertaken  to  explore  the  effects  of  alloy  formation  on 
the  surface  properties  of  the  binary  constituents.  A  portion  of  this  data 
will  be  presented  here  to  illustrate  central  aspects  of  the  oxidation,  such 
as  the  observed  instability  of  the  oxide  formed;  a  more  complete  exposition 
of  the  findings  will  be  published  seperately. 

II.  Experimental 

Samples  of  HgTe  and  Hgg.ggCdr). 31  Te ,  grown  by  solid  state  recrystalli¬ 
zation,  were  transferred  through  an  air  lock  into  the  vacuum  system  after 
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baking  to  avoid  heating  the  specimens.  Clean  surfaces  of  these  materials 
and  a  bar  of  CdTe  (II-VI,  Inc.)  were  prepared  by  cleaving  in  vacuum  (p  ~ 

5  x  10-11  torr)  prior  to  initiating  the  oxidation  sequence. 

The  oxidation  was  carried  out  in  the  small  evacuated  volume  from 
which  the  samples  had  been  originally  transferred  after  separately  baking 
this  portion  of  the  apparatus  to  reduce  the  pressure  and  water  content. 
In  this  respect  the  oxidation  process  was  similar  to  that  reported  previ¬ 
ously  [7],  To  stimulate  the  otherwise  extremely  slow  uptake  of  oxygen 
[7]  an  unbiased  hot  thoriated  iridium  filament  was  positioned  in  line-of- 
sight  of  the  sample  surface  during  the  exposure.  The  oxygen  pressure  was 
monitored  by  a  millitor  ionization  gauge  (Varian  Associates)  not  in 
direct  1 i ne-of-sight  of  the  sample,  and  the  pressure  regulated  at  a  fixed 
value  of  5  x  10“^  Torr.  Before  beginning  the  oxidation  study,  it  was 
ascertained  that  only  the  combination  of  oxygen  and  the  hot  filament  in 
front  of  the  sample  resulted  in  the  reactions  reported;  operation  of  the 
out-of-1  i ne-of-si ght  gauge  in  oxygen  or  use  of  the  hot  filament  in  the 
absence  of  an  oxygen  ambient  produced  no  detectable  change  in  the  state 
of  the  surface. 

The  oxidation  procedure  in  the  present  study  is  distinguished  from 
an  earlier  one  [7]  in  the  establishment  of  a  linear  exposure  scale.  This 
was  accomplished  by  repeatedly  applying  a  fixed  dose  (5  x  10"5  Torr  for  500 
sec.)  rather  than  manipulating  the  pressure  to  increase  the  exposure.  In 
the  latter  procedure,  the  pressure  variation  of  the  surface  reactions 
reflects  both  the  dependence  of  the  oxidation  on  exposure  as  well  as  the 
pressure  dependent  efficiency  of  the  uptake  activation  process.  By  repe¬ 
tition  of  a  fixed  dose  instead,  a  sample  expi  sed  twice  to  the  activated 
oxygen  received  double  the  dose  as  one  exposed  once  only.  The  linear 


90 


} 

RD-A156  079  HGCDTE  SURFACE  AND  DEFECT  5TUDV  PROGRAM(U)  SANTA 

BARBARA  RESEARCH  CENTER  GOLETA  CALIF  J  A  WILSON  ET  AL. 
01  JUL  84  RPT50220  MDA902-82-C-0108 
UNCLASSIFIED  F/G  28/12 


exposure  scale  established  in  this  way  facilitates  preliminary  analysis 
of  the  kinetics  involved  in  the  oxygen  uptake. 

Photoemission  spectroscopy  served  to  characterize  the  clean  surface 
and  monitor  the  effect's  of  oxidation.  To  excite  the  photoelectrons, 
photon  energies  of  55  and  85  eV  from  a  synchrotron  source  (Stanford 
Synchrotron  Radiation  Laboratory)  were  used  for  spectra  of  the  Hg  5d  and 
Cd  4d  core  levels  and  valence  band  and  the  Te  4d  core  lines,  respec¬ 
tively.  This  choice  provides  high  surface  sensitivity  and  results  in 
nearly  the  same  kinetic  energy  and  thus  inelastic  mean  free  path  for  the 
core  level  signals  of  all  substrate  constituents.  The  monochromator  was 
operated  with  0.2  eV  resolution  and  the  double  pass  cylindrical  mirror 
analyzer  used  to  energy  analyze  and  collect  the  electron  distribution 
gave  an  energy  resolution  of  0.18  eV.  The  spectra  presented  have  been 
normalized  to  the  incident  beam  intensity  as  monitored  by  the  electron 
yield  from  an  80  percent  transmitting  copper-coated  stainless  steel  mesh 
inserted  into  the  beam. 

III.  Results  and  discussion 

The  chemistry  and  kinetics  of  the  oxidation  reactions  are  revealed 
in  shifts  in  energy  and  changes  in  intensity  of  the  core  line  emission  as 
a  function  of  activated  oxygen  exposure.  The  Cd  4d  and  Hg  5d  core  lines 
and  valence  band  emission  from  the  clean  and  oxygen  exposed  surface  of 
Hgo.69^^0. 31  Te ,  shown  in  Figure  1,  illustrate  key  features  of  the  alloy 
surface  oxygen  chemistry.  As  the  number  of  doses  increases,  an  oxygen  2p 
derived  structure  appears  and  systematical  ly  increases  in  strength  near 
46  eV,  indicating  the  uptake  of  oxygen  by  the  surface.  As  the  oxygen 
coverage  increases,  the  spin  orbit  split  Hg  5d 5/2  ancl  ^3/2  levels  are 
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attenuated  while  the  Cd  4d  emission  is  substantially  unchanged  in  intensity. 
The  difference  in  the  exposure  dependence  of  the  Hg  and  Cd  signals  indicate 
the  oxide  grows  excluding  Hg  while  incorporating  Cd.  The  release  of  Hg 
from  the  growing  film  is  similarly  found  to  occur  upon  exposure  of  HgTe 
to  activated  oxygen  [8].  In  the  limited  number  of  exposures  performed  on 
CdTe,  Cd  is  observed  as  a  constituent  of  the  oxide  formed  on  this  compound. 

The  Te  4d  core  level  spectra  as  a  function  of  exposure,  presented  in 
Figure  2,  augment  the  information  gained  from  the  cation  core  lines  in 
describing  the  oxidation  process.  Characteristic  of  oxidized  Te  [9],  a 
doublet  shifted  by  3.6  eV  to  higher  energy  appears  with  the  first  dose  and 
increases  in  strength  with  subsequent  exposure  (Figure  2).  This  chemical 
shift,  the  difference  in  energy  between  the  peaks  of  the  oxidized  Te  4d5/£ 
signature  and  the  corresponding  emission  from  Te  in  the  semiconductor 
lattice,  is  initially  the  same  for  HgTe  and  the  alloy.  Additional  exposure 
results  in  a  different  chemical  shift  for  the  two  materials  and  produces 
an  oxidized  Te  signal  broadened  by  0.2  eV  (full  width  at  half  maximum) 
for  both  samples  (see  Figure  2  for  the  alloy  case.)  These  changes  in  peak 
position  and  width  of  the  oxidized  Te  reflect  a  changing  mixture  of 
oxides,  each  component  of  which  has  an  oxidized  Te  signal  at  somewhat 
different  energy. 

That  the  oxidized  Te  signal  contains  more  than  one  contribution  is 
more  directly  illustrated  in  Figure  3a,  which  also  indicates  the  temporal 
stability  of  these  compounds.  The  two  spectra  characterize  the  HgTe 
surface  after  the  highest  exposure  used  and  were  recorded  immediately 
following  the  dose  and  12  hours  later.  With  the  passage  of  time,  the  4d 
peaks  arising  from  Te  in  the  oxide  sharpen  and  appear  at  lower  binding 
energy  relative  to  the  core  level  emission  immediately  following  the 


exposure.  It  thus  appears  that  the  oxide  formed  at  the  higher  exposure 
levels  includes  a  component  (indicated  by  the  higher  binding  energy 
contribution)  which  is  unstable  at  room  temperature.  Interfacial  reac¬ 
tions  which  convert  one  type  of  oxide  into  another  have  been  predicted 
from  thermodynamic  data  and  observed  between  bulk  compounds  under  moder¬ 
ate  heating  [10];  the  oxides  contributing  to  measured  spectra  (see 
Figure  3)  must  be  identified  before  the  predicted  reactions  can  be  veri¬ 
fied  in  this  case.  The  sharpening  shown  in  Figure  3a  has  been  observed 
in  a  previous  study  of  the  growth  of  thin  native  oxide  layers  on  HgCdTe 
[7]  but  could  not  be  assigned  unambiguously  to  an  oxide  instability  at 
that  time.  A  similar  narrowing  and  sharpening  of  the  oxide  signal  was 
observed  when  the  HgCdTe  sample  was  re-examined  23  hours  after  completing 
the  highest  dose.  This  phenomena  was  not  examined  on  CdTe. 

Like  the  oxidized  Te  linesnape  discussed  above,  the  Te  4d  emission 
near  40.8  eV  in  Figure  2  from  the  oxygen  exposed  surface  also  arises  from 
more  than  one  contribution.  In  the  spectrum  taken  following  the  first 
exposure,  the  Te  4d  signal  at  the  bulk  (clean  surface)  position  is  broad¬ 
ened  asymmetrically  to  higher  binding  energy.  This  broadening  is  also 
apparent  following  additional  exposure,  particularly  after  doses  3  and  4, 
but  tends  to  decrease  overall  with  continued  uptake.  As  illustrated  in 
Figure  3b,  the  additional  emission  at  higher  binding  energy  can  be 
accounted  for  by  a  replica  of  the  clean  surface  peaks  shifted  by  0.6  eV 
to  higher  binding  energy  and  added  to  the  emission  from  the  Te  in  the 
HgCdTe.  The  new  state  of  Te  has  a  binding  energy  characteri stic  of 
elemental  Te  [9]  although  a  partially  oxidized  Te  signal  with  the  same 
binding  energy  would  not  be  distinguished  from  elemental  Te  in  these 
measurements.  The  Te  signal  at  the  elemental  position  was  found  in 


varying  amounts  on  ail  three  compounds. 

In  addition  to  the  chemical  information  derived  from  the  photo 
emission  spectra,  preliminary  kinetic  data  can  also  be  extracted. 

Assuming  an  abrupt  interface  between  the  oxide  film  and  the  underlying 
semiconductor,  an  equal  density  of  Te  atoms  in  the  two  regions,  and  a 
constant  inelastic  mean  free  path  in  each  case,  the  thickness  of  the 
overlayer  may  be  estimated  from  the  exponential  attenuation  of  the  sub¬ 
strate  signal  with  oxide  thickness  and  the  increasing  strength  of  the 
emission  from  oxidized  Te  (Figure  2).  Working  in  this  simplified  model, 
the  thickness,  in  units  of  the  inelastic  mean  free  path,  X,  is  given  by 


where  [Te0x]  and  [Te]  are  the  integrated  intensities  of  the  core  level 
emission  from  the  oxidized  and  unoxidized  Te  (including  the  signal  at  the 
elemental  position).  The  results  of  such  an  analysis  for  the  three 
samples  studied  is  presented  as  a  function  of  the  number  of  standard 
doses  applied  in  Figure  4.  Based  on  this  measure  of  thickness,  the 
oxidation  is  found  to  proceed  in  proportion  to  the  exposure  in  each  case 
with  no  evidence  in  the  data  for  saturation  in  the  uptake  for  films  grown 
in  excess  of  at  least  two  monolayers  thick.  The  rate  of  oxidation,  given 
by  the  slope  of  the  lines,  is  found  to  be  69  percent  as  great  for  the 
h90.69^0. 31 Te  as  for  the  HgTe  sample;  i.e.,  the  rate  of  oxidation 

scales  as  the  number  of  Hg-Te  bonds.  This  finding  in  conjunction  with 
the  result  from  the  core  level  spectra  that  Hg  is  released  and  not  incor¬ 
porated  in  the  initial  stages  of  oxide  growth  suggests  the  importance  of 
breaking  the  Hg-Te  bond  as  a  step  in  the  oxidation  process.  If  HgCdTe 
were  to  oxidize  as  merely  a  mixture  of  the  binary  compounds,  the  rate  for 
for  the  alloy  would  be  an  appropriate  average  of  the  rates  for  HgTe  and 


CdTe  and  the  line  for  the  alloy  in  Figure  4  would  fall  intermediate  to 
those  of  the  binaries.  Instead,  the  rate  of  oxidation  as  determined  from 
the  change  in  thickness  with  exposure  is  greater  for  CdTe  than  the  x  =  0.31 
alloy,  implying  a  modification  of  the  activity  of  the  Cd-Te  site  in 
^90.69^0. 31  Te  compared  to  the  binary. 

IV.  Conclusion 

The  spectra  obtained  during  study  of  the  activated  oxidation  of  CdTe, 
HgTe,  and  Hgo^sgCdo.si  Te  indicate  the  oxide  which  initially  forms  on  the 
samples  containing  Hg  does  so  with  the  release  of  this  constituent; 
whereas,  Cd  is  retained  in  the  oxide  of  CdTe  and  the  alloy.  The  observa¬ 
tion  of  the  release  of  Hg  together  with  the  rate  limiting  role  of  the 
Hg-Te  bond  are  consistent  with  the  lesser  strength  of  this  bond  compared 
to  that  of  Cd  to  Te.  Similarly,  the  modification  in  the  activity  of  the 
Cd-Te  site  in  the  alloy  exemplified  by  the  higher  rate  of  oxidation  for 
CdTe  than  the  x  =  0.31  sample  is  in  line  with  the  difference  in  elec¬ 
tronic  structure  of  the  two  materials  which  results  in  a  stronger  Cd-Te 
bond  in  the  alloy  than  CdTe.  It  would  be  of  interest  to  extend  the 
measurements  to  higher  x  value,  where  the  kinetics  operative  during  the 
oxidation  of  CdTe  may  become  dominant.  Further  study  of  the  kinetics  of 
oxidation  is  also  warranted  in  as  much  as  the  oxide  produced  by  the 
activated  oxidation  procedure  on  HgTe  and  HgCdTe  exhibited  unstable 
behavior  indicative  of  nonequilibrium  constitution  or  interfacial  reac¬ 
tions. 
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Figure  1 


Figure  2 
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Figure  4 
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Photoemi ssion  spectra  of  the  Cd  4d,  Hg  5d  and  valence  band 
emission  of  HgQ#5gCdo,3i Te  cleaved  in  vacuum  and  exposed 
repeatedly  to  fixed  doses  of  oxygen  activated  by  a  hot 
filament  in  1 ine-of-sight  of  the  sample  surface. 
Photoemission  spectra  of  the  Te  4d  lines  on  the  clean  and 
oxygen  dosed  surface.  The  details  of  the  changes  noted  in 
the  oxidized  Te  and  unoxidized  signal  are  discussed  in  the 
text. 

Illustration  of  time  evolution  at  room  temperature  of  the 
oxidized  Te  4d  signal  following  the  highest  exposure  on 
HgTe.  After  12  hours,  the  peaks  sharpen,  losing  the  high 
binding  energy  component.  A  similar  result  occurred  for 
the  H90.69Cd0.31  Te- 

Decomposition  of  the  Te  4d  emission  from  HgQ^ggCdf)^ 31  Te 
recorded  following  the  first  exposure  into  peaks  at  the 
energy  of  Te  in  HgCdTe  and  a  contribution  at  the  energy 
of  elemental  Te. 

Thickness  of  the  oxide  derived  from  the  relative  intensity 
of  the  oxidized  and  unoxidized  Te  emission  as  a  function  of 
the  number  of  standard  doses  applied  to  the  surfaces  of 
CdTe,  HgTe,  and  Hgg.ggCdQ^ 31 Te.  One  escape  length,  the 
unit  of  thickness,  is  equivalent  to  one  to  three  mono- 
1 ayers. 
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d  =  (l  -x)do(AC)  +  xd0(BC) 


(9) 
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where  d^AC)  and  d^BC)  are  the  pure  crystal  lattice  constants.  Then  in 
an  A(3)B(l)  configuration  [see  Fig.  1(a)],  the  anion  will  be  displaced  along 
the  BC  bond.  For  definiteness,  suppose  d(AC)  >  d(BC);  then  the  anion  will 
be  shifted  away  from  the  tetrahedral  site  toward  the  B  atom.  Because 
d^AC)  s  3  5r  d0(BC),  the  d(BC)  bond  will  be  elongated  and  the  d(AC)  bonds 
will  also  be  locally  elongated  in  a  strained  configuration.  We  can  for  this 
case  [see  Fig.  1(a)]  write 


d(BC)  =  (1  -  <5)d 


(10) 


and  for  the  three  AC  bonds 


d(AC)  =  V1  +  Tp+<52  3  (ID 

In  an  A(2)B(2)  configuration  [see  Fig.  1(b)],  we  have  for  both  BC  bonds 

d(BC)  =  U2  d  (12; 


and  for  both  AC  bonds 


d(AC)  -  V 1  +  r5  +  a 

T:''h  b’-  ‘.he  •  onddion  dF-/  de  =  0 


3 


The  values  of  6  r 


or 


!  _  StU^-Uj-tBC)^-]2  _  2[ug-Uj-(AC)v1-]; 


AEb1  = 


tr§  -  ei(BC) 


eb°  -  ei(AC) 


dd  2  (  [Ug-Uj-tBOV!,]2  _  [tlg-U^AC)^.]2 
4dd  *  <  e£  -  £b  (BC)  eb°  -  c.1  (AC) 


(?) 


[UbL  (BC)  U^Vi-f  [UbLCACjU^-V!-]2  ) 
cb  (BC)  -  e°  Eb  (AC)  -  £a°  ' 

where  the  partial  derivative  is  taken  with  respect  to  the  central  bond 
length  (a  term  designated  with  a  superscript  o)  and  in  this  simple 
approximation 


AE-?  =  tjE-J  ,77  0, 1,2,3  .  (8) 

The  bond  length  shifts  are  found  from  Eq.  (2).  The  modification  of  these 
expressions  for  an  anion-substituted  alloy  follows  from  symmetry. 

Method.  2— A  more  realistic  calculation  includes  strain  and  permits 
the  lattice  to  fit  together.  The  simplest  non-self-consistent,  near¬ 
neighbor  version  of  this  calculation  proceeds  as  follows.  In  a  cation- 
substituted  alloy  \;_XBXC  (e  g.  Hgi_xCdxTe  or  Zn]_xHgxTe),  the  cations  are 
assumed  to  occupy  their  regular  fee  subiattice  sites,  while  the  amor;.' 
accommodate  to  .=  ,r  local  configuration  The  average  lattice  cons  lam 
taken  to  follow  s  rule 
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fits  together  by  having  the  metal  sublattice  nearly  retain  its  f.c.c.  struc¬ 
ture  with  a  lattice  constant  following  Vegard’s  rule,  and  each  As  atom 
adjust  its  local  position  and  accommodates  to  which  metal  atoms  occupy 
its  four  neighboring  sites.  Thus  there  are  five  local  arrangements;  in  the 
first  an  As  atom  is  surrounded  by  four  Ga  atoms  [a  Ga(4)In(o) 
configuration],  in  the  second  there  are  three  Ga  atoms  and  one  In  [a 
Ga(3)In(l)  configuration],  and  the  last  three  configurations  are 
Ga(2)In(2),  G(l)In(3),  and  Ga(o)In(4).  In  a  general  cation  Aj-^B^C  alloy  or 
anion  CAi_xB*  alloy,  the  generalizations  are  obvious,  and  one  can  examine 
A(4  —  Tj)B(rj);  77  =  1,2, 3, 4  configurations. 

Method.  1— The  lowest  level  of  approximation  that  takes  account  of 

these  findings  calculates  [from  Eqs.  (1)  and  (2)]  a  given  bond's  energy 

« 

and  length  shift  when  its  surrounding  bonds  have  different 
constituents. The  calculation  is  done  assuming  that  the  surrounding  bonds 
retain  their  pure  crystal  bond  lengths.  This  model  is  somewhat  unrealis¬ 
tic,  in  that  it  ignores  the  strains  that  must  be  present  to  allow  the  lattice 
to  fit  together.  However,  it  is  nonetheless  useful  because  it  treats  the 
largest  terms,  viz.  those  arising  from  charge  transfer  between  the  substi¬ 
tuted  species  in  the  alloy,  and  results  in  analytic  expressions,  while  the 
successively  more  complicated  treatments  must  be  done  numerically.  It 
is  also  useful  later  because,  taken  together  with  the  more  complete  cal¬ 
culations,  it  allows  one  to  judge  the  relative  contribution  from  different 
physical  effects  From  Eq.  (1)  in  a  cation-substituted  alloy  one  can 
deduce  an  expression  for  the  energy  shift  of  an  AC  bond  when  p  sur¬ 
rounding  bonds  are  replaced  by  BC,  AE|  as 
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Finally  combining  Eqs.  (1)  through  (5),  the  bond  energy  is 


Eb  =  -(1  +  a*)(Vf  +  Vf)1/2  -  |-(1  -  a2)2 


(v?+  +  VfJ 

(V|  +  V|)1/2 


+  (V1+  +  V,_i  +  hUZJsz.  .  (6) 

Equation  (6)  was  first  derived  by  Harrison  using  a  much  simpler  argu¬ 
ment.1  The  bond  energies  quoted  in  Table  I  and  on  Fig.  5  are  calculated 
from  this  expression,  using  as  input  the  correlated  term  values  from 
Table  I  and  the  experimental  bond  lengths. 

B.  Alloy  Calculations 

There  are  several  layers  of  sophistication  that  can  be  used  to  com¬ 
pute  the  bond  energy  and  length  changes  between  pure  compounds  and 
their  alloys.  However,  all  the  calculations  are  motivated  by  the  EXAFF 
experimental  result  found  first  for  Gai_xInxAs  and  later  in  other  com¬ 
pounds.5  It  was  found  that,  while  the  average  bond  length  in  the  alloy  fol¬ 
lows  Vegard’s  rule  and  varies  Linearly  between  those  of  the  pure  consti¬ 
tuents,  the  individual  bond  lengths  of  GaAs  and  InAs  are  nearly  the  same 
as  those  of  the  pure  materials—changing  by  only  about  4  percent.  The 
shorter  GaAs  bond  increases  its  length  when  it  is  an  impurity  in  an  InAs 
host  by  about  4  percent,  and  the  longer  InAs  bond  shrinks  by  about  4  per¬ 
cent  in  the  opposite  extreme  where  it  is  the  impurity  in  GaAs  The  lattice 


(4) 


\ 

s  £  (Tm+.n  4*  Tjn-n) 
i«0 


The  last  identity  in  Eq.  (4)  defines  the  indicated  bond  tensions.  Collecting 
all  the  results  from  the  earlier  paragraphs,  one  finds  (dropping  the  bond 
designation  n): 
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Harrison  has  written  Eq.  (1)  without  explicit  reference  to  the  sur¬ 
rounding  bonds.  Equation  (l)  in  its  present  form  is  rather  easily  used  for 
alloy  calculations  by  replacing  host  atom  terms  by  one  appropriate  to  the 
alloy. 

The  parameter  y  for  each  compound  is  fixed  by  the  condition  that 
there  is  no  net  force  on  atoms  at  their  equilibrium  lattice  positions.  This 
condition  is  satisfied  if  the  gradient  of  the  total  energy  with  respect  to  all 
the  bond  lengths  vanishes: 


Vfdij  Et  =  0,  where  ET  =  £  E£ 

n 


This  condition  is  equivalent  to  the  requirement  that  the  net  bond  tension 
Ton  +  Trm  for  the  nth  bond  equals  the  repulsive  force 


T  +  T  —  « 

1on  lmn  ,5  ~  u 

Qn 


(2) 


where 


Ton  =  J“[— 2- Vvf  +  Vf  j  =  4ac  Vvf  +  V$/ d 
ddn 


(3) 


T„ 


=  2 


dd- 
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■s 

:=0 


rr.+  .n 


^  T. 
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lil 


■p 


are  the  probability  amplitudes  of  finding  the  cation  and  anion  sp3  -hybrid 
contributions  to  the  bonding  and  antibonding  states,  e.g.  Ub+  is  the  pro¬ 
bability  amplitude  of  finding  the  cation3  hybrid  in  the  bonding  state,  etc., 

ap  =  V3/  W|  +  Vf  is  the  polarity, 


ac  s  |  V2  |  /  Vvf  +  vf  is  the  covalency, 


«,b,  =  1*1  VvJTv? 


are  the  bonding  and  antibonding  one  electron  energy  states  owing  only  to 
covalent  and  ionic  interactions, 


€  rr. ' 1 1 


i  -  1,2,3 


are  the  contribir  -  ‘o  the  bond  in  question  arising  from  its  antibond;::^ 
state  interaction.-  •  th  a  neighboring  bonding  state  1,  and  y  is  an  adju¬ 
stable  coefficient  >■':  »o:i  ui  the  measured  bond  lengths  are  reproduced 
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transfer  from  the  cation  to  the  anion),  the  second  term  is  the  promotion 
energy  per  bond  to  form  sp3  hybrids  on  both  atoms,  the  third  term  is  the 
bond  formation  energy  owing  to  the  covalent  and  ionic  terms,  the  fourth 
term  is  Harrison's  metallization  energy,  and  the  fifth  term  is  the  repul¬ 
sive  overlap  energy  that  prevents  the  bonds  from  collapsing.  The  various 
symbols  are  defined  below: 

cp|±j  is  the  {  }  p-state  correlated  term  values4 

V„±j  s  (ep(±|  -  eS|±j )  /  4  are  hopping  integrals  between  two  adjacent 
bonds  coupled  through  a  l  aniQn  j 

c9(±j  are  the  {  $  s-state  correlated  term  values 

V%  s  -  24.5 /d2  [eV]  is  the  covalent  energy 
d  is  the  bond  length 

Vas(e+-0/2  IS  P°,<ar  energy 
«|±|  =  (e3|±j  +  3epj±])/4  are  the  sp3  hybrid  energies 

c&W  ■  -  ei) 

j=i 


is  the  contribution  to  the  reduction  of  the  bonding  state  energy  e°  calcu¬ 
lated  by  second-order  perturbation  theory  of  the  bond  in  question 
(denoted  by  superscript  o),  owing  to  its  interaction  with  the  three  neigh¬ 
boring  (i  =  1.2,3)  antibonding  states  sharing  the  same  [  J  } 
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expected,  because  the  electron  effective  mass  in  a  narrow-gap  material  is 
small;  hence,  the  electron  wave  functions  at  Tic  are  distributed  over 
many  atoms,  with  the  result  that  alloy  scattering  rates  are  small  at  the 
band  edge.  The  hole  mobilities  are  large  for  a  somewhat  different  rea¬ 
son:  The  valence  band  edge  is  mostly  composed  of  Te  p-  states;  however, 
there  is  also  some  cation  p-  state  contribution.  Because  the  p-  state 
energies  of  Hg  and  Zn  are  nearly  the  same,  the  alloy  scattering  strength 
at  the  valence  band  edge  is  small.  We  conclude  that,  because  of  its  struc¬ 
tural  advantage,  HgZnTe  may  be  superior  to  HgCdTe  for  infrared  devices. 


II  ALLOY  MODIFICATION  OF  THE  BOND  ENERGY  AND  LENGTH 

A.  Pure  Compounds 

The  first  task  is  to  recast  Harrison’s  bond  energy  and  bond  length 
formalism  into  a  structure  suitable  for  generalization  to  alloys  Focus 
attention  on  the  nth  bond  of  a  pure  zinc  blende  structured  compound 
semiconductor;  for  the  time  being,  suppress  any  notation  identifying  it 
Then,  in  a  II-VI  compound,  the  bond  energy  is  is:1 

Eb  =  Cp*  ~  C?~  +  VI+  +  Vj_  !  — 2(Vf  +  V3*)  1/2  +  2  V  (e^.  +  . 

2  i=o  d 

where  the  first  term  is  the  energy  per  bond  needed  to  transfer  two  elec¬ 
trons  from  a  p-  .stale  on  the  Group  VI  atom  (the  anion)  to  a  p-  state  on 
Mm  the  Group  1!  atom  (the  cation),  so  that  both  start  with  four  elec¬ 
trons  (ultimately  in  the  final  bonding  arrangement  there  is  a  net.  electro:: 
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I  INTRODUCTION 

The  objective  of  this  work  is  to  understand  the  microscopic  mechan¬ 
isms  that  govern  the  stability  of  Hgi_xCdxTe  alloys,  and  then  to  suggest 
changes  in  the  material  that  improve  its  strength  without  adversely 
affecting  its  electronic  behavior.  In  pursuit  of  this  goal  we  have  extended 
Harrison’s  bonding  theory,1  which  is  applicable  to  all  tetrahedral  struc¬ 
tured  semiconductors,  to  calculate  bond  length  and  energy  changes  in  an 
alloy— including  charge  shift  and  reconstruction  effects— relative  to  their 
pure  crystal  values,  vacancy  formation  energies,  dislocation  energies, 
and  hardness. 

In  this  paper,  we  review  the  experimental  situation  in  several  of 
these  areas  and  compare  some  of  the  results  with  theory  for  all  the 
Group  IV,  III-V  compound,  and  II-VI  compound  semiconductors  and  their 
alloys.  For  example,  the  theory  properly  predicts  the  observed  inverse 
ninth  power  bond  length  (d-9)  dependence  of  the  hardness2  of  semicon¬ 
ductors  and,  with  no  adjustable  parameters,  their  correct  magnitude.  We 
demonstrate  that  the  weak  HgTe  bond  is  destabilized  by  alloying  it  with 
CdTe,  HgS,  or  HgSe;  however,  the  bond  is  stabilized  by  ZnTe.  Moreover, 
because  the  bond  length  of  ZnTe  (2.406  A)  is  14  percent  shorter  than  that 
of  HgTe  (2.797  A)  or  CdTe  (2.804  A),  the  dislocation  energy  per  unit  length 
and  hardness  of  the  alloys  Hg2_xZnxTe  and  Cdj_xZnxTe  are  predicted  in 
agreement  with  some  experiments  to  be  significantly  higher  than  those  of 
the  compounds  with,  x  =  0  Measurements  indicate  that  the  electron  and 
hole  mobilities  of  !>g._xZnxTe  with  x«.0.16,  corresponding  to  a  0.1-eY 
band  gap,  are  comparable  to  those  of  Hgj_xCdxTe  for  x  =  0  2. 3  This  is  to  be 
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improve  the  structural  properties  of  HgTe  and  CdTe.  Experiments  that 
support  these  predictions  are  noted.  The  electronic  transport  properties 
of  0.1-eV  band  gap  HgZnTe  are  about  the  same  as  those  of  HgCdTe,  and 
the  structural  properties  of  the  Zn  compound  are  superior;  thus,  we  con¬ 
clude  that  HgZnTe  is  likely  to  be  the  better  material  for  IR  devices. 


Section  IV-I 
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(14) 


Equation  (14)  is  solved  by: 

(1)  Assigning  the  -yns  their  pure  crystal  values 

(2)  Supposing  that  the  four  tetrahedrally  configured  anions  are 
embedded  in  the  virtual  crystal  [see  Fig.  1(c)] 

(3)  For  each  x,  <5  is  varied  and  the  bond  tensions  are  calculated  and 
inserted  into  Eq.  (14)  until  a  value  of  6  is  found  that  satisfies 
this  equation. 

The  computer  outputs  are  <5  and  the  AC  and  BC  bond  energies.  This 
method  can  be  extended  to  take  explicit  account  of  the  bond  length  vari¬ 
ations  of  larger  clusters,  and  3  can  be  calculated  explicitly  (instead  of 
using  Vegard’s  rule)  so  the  results  can  be  made  to  be  self-consistent 
Also,  effects  arising  from  bond  angle  distortions  can  be  included.  The 
major  effects  that  are  currently  included  in  the  calculation  arise  from 
the  bond  length  dependence  of  V2,  and  charge  shifts  driven  by  hybrid 
energy  differences  and  coupled  through  the  metallization  terms.  How¬ 
ever,  V2  is  also  sensitive  to  bond  angle  distortions,  and  these  effects  are 
not  included  in  the  present  results.  We  expect  their  inclusion  will  modify 
the  quantitative  results  by  10  to  30  percent,  but  not  the  trends. 


C.  Results 

The  correlated  atomic  term  values4  for  the  elements  in  the  more 
common  II- VI  compounds  are  listed  in  Table  I,  along  with  the  bond  ener¬ 
gies  of  their  nine  >  umpounds.  The  valence  s  -state  energies  for  Zn  are 
large,  become  smaller  as  for  Cd,  and  then  (in  an  unusual  occurrence)  the 
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trend  reverses  and  the  Hg  s-  levels  are  deeper  again.  This  is  caused  by 
the  relativistic  terms  that  become  important  in  Hg.  The  p-  state  levels 
have  a  similar  trend;  however,  the  differences  among  the  atoms  are  much 
smaller.  The  levels  of  the  anions  are  deeper  than  those  of  the  cations. 
Thus,  there  is  a  net  electron  transfer  from  the  cations  to  the  anions, 
which  is  responsible  for  the  ionic  contribution  to  the  bonding.  The  polar¬ 
ity,  also  listed  in  Table  I,  reflects  the  relative  contribution  of  ionic  and 
covalent  character  to  the  bonds.  Notice  that  ap  is  small  for  ZnTe  larger 
for  CdTe,  and  smaller  again  for  HgTe.  By  contrast,  the  metallic  contribu¬ 
tion  Em  is  large  for  ZnTe  and  HgTe  and  small  for  CdTe.  In  fact,  HgTe 
would  not  be  bound  if  it  were  not  for  the  metallic  contribution  to  the 
bond,  not  a  terribly  surprising  result,  because  HgTe  is  a  semimetal.  The 
bond  lengths  of  CdTe  and  HgTe  are  very  nearly  the  same,  but  this  is  an 
accident  resulting  from  the  balance  of  contributions  to  their  bonds: 
CdTe  is  dominated  by  the  covalent  and  ionic  contributions  to  its  bond, 
while  HgTe  has  bonds  that  are  mostly  covalent  and  metallic.  We  should 
therefore  not  be  too  surprised  if  their  alloy  displays  some  unusual  bond 
length  variations. 

Figures  2  and  3  illustrate  the  Hgj_xCdxTe  alloy  concentration  varia¬ 
tions  of  the  length  and  energy  for  the  five  configurations.  These  results 
were  calculated  using  the  second  method  discussed  in  Section  II-B.  In 
Fig  2  we  see  the  unexpected  result  that  the  shorter  HgTe  bond  becomes 
still  shorter  in  the  alloy,  while  the  CdTe  bond  length  increases.  Moreover, 
the  changes  are  large  compared  to  the  original  length  difference  Thus 
one  can  expect  to  find  local  microstrains,  even  in  this  nearly  1  a  1 1  ■  -  -  ■  - 
matched  material  The  bond  energy  variations  shown  in  Fig  3  ;  - 
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display  an  important  result.  The  already  weak  HgTe  bond  is  destabilized 
by  the  presence  of  Cd.  The  flat  concentration  variation  of  the  bond 
lengths  and  energies  is  peculiar  to  HgCdTe;  much  more  structured 
behaviors  are  exhibited  by  other  compounds. 

Several  considerations  to  keep  in  mind  when  dealing  with  alloys  are 
depicted  schematically  in  Fig.  4  for  a  Hg(3)Cd(i)  configuration.  Due  to 
the  charge  shift  alone,  the  HgTe  bond  shrinks  by  0.030  A,  the  minimum  in 
the  configuration  diagram  moves  to  a  smaller  bond  length  and  the  depth 
of  the  well  is  decreased.  Then  the  bond  is  stretched  by  0.011  A  to  its  final 
length.  Thus,  the  bond  is  not  centered  at  a  local  minimum.  For  example, 
this  implies  that  the  force  constant  of  local  phonon  modes  will  be  related 
to  the  second  derivative  of  this  configuration  diagram  evaluated  at  the 
strained  position  away  from  the  well  minimum.  This  effect  is  even  more 
pronounced  for  materials  with  a  larger  lattice  mismatch.  Notice  that  the 
CdTe  bond  is  also  stretched  relative  to  its  local  minimum  in  this 
configuration. 

The  result  for  bond  energy  variations  for  all  the  II-VT  compound 
alloys  is  gathered  in  Fig.  5.  Here  the  energy  shifts  of  an  AC  bond  in  a 
A(l)B(3)  configuration  AE£  =  3AE.J  are  presented  as  a  percentage  of  the 
pure  material  bond  energy  for  cation  substitute  alloys  along  with 
equivalent  results  for  anion  substantial  alloys,  e  g  HgSj.jTe*.  The  energy 
shifts  are  calculated  following  the  first  method  presented  in  Section  II-B 
We  have  not  had  the  opportunity  to  run  all  the  compounds  following  the 
more  sophisticated  calculational  procedure  The  bond  at  the  top  of  a 
column  is  the  one  whose  energy  shift  is  calculated,  and  the  other 


constituent  is  identified  in  the  line  that  indicates  the  energy  shift.  A 
negative  3AEb/E},  is  destabilizing  and  a  positive  one  stabilizing.  Focus  on 
HgTe.  It  is  destabilized  by  CdTe,  HgS,  and  HgSe,  but  is  stablized  by  ZnTe. 
The  CdTe  destabilizes  the  HgTe  bond  by  18  percent  according  to  the  sim¬ 
ple  calculation  and  by  15  percent  according  to  the  better  one.  The 
difference  between  Zn  and  Cd  occurs  because  in  the  Cd  case  there  is  a 
net  electron  transfer  from  the  Cd  to  the  Hg,  while  in  the  case  of  Zn  there 
is  practically  no  such  transfer.  These  extra  transferred  electrons  must 
occupy  antibonding  states,  because  the  bonding  states  on  the  HgTe  are 
full,  so  the  net  effect  (including  other  energy  state  shifts)  is  to  destabilize 
the  HgTe  bond.  This  predicted  destabilization  of  the  HgTe  bond  in 
HgCdTe  alloys  agrees  with  the  d-  state  energy  shift  reported  by  Spicer  et 
al.  at  last  year’s  workshop,®  and  by  the  oxygen  uptake  experiments 
presented  by  the  same  group  at  this  conference.7 


IQ  DISLOCATION  ENERGIES  AND  HARDNESS 

The  dislocation  energies  and  hardness  of  semiconductor  are  calcu¬ 
lated  by  an  extension  of  Harrison's  method.1  It  is  demonstrated,  in  agree¬ 
ment  with  experiment,  that  dislocation  energies  per  unit  length,  depend¬ 
ing  on  the  covalency,  are  proportional  to  d-9,  where  d  is  the  bond  length 
and  hardness  is  proportional  to  d-5  to  d-n  The  low  powers  are  for  pure 
covalent  material.'  The  hardness  is  related  to  the  interaction  energic - 
among  dislocattor  -  The  detailed  theory  will  be  published  elsewhere  8 


Figure  4  shows  the  theoretical  hardness  of  several  semiconductors 
plotted  against  experimental  values.8  The  theory  is  calculated  for  two 
different  dislocation  configurations  caused  by  the  extra  material  pushed 
by  the  indentor  into  the  crystal.  The  lower  end  of  each  vertical  line  is  a 
realistic  hardness  for  the  minimum  energy  configuration,  and  the  upper 
ends  of  the  lines  are  the  values  found  for  a  higher-energy  configuration. 
In  most  cases,  one  would  expect  real  configurations  to  have  hardnesses 
lying  somewhere  between  these  two  values,  once  the  effect  of  heat  gen¬ 
erated  by  the  indention  process  is  subtracted  from  the  experimental 
numbers.  This  theory  has  no  adjustable  parameters,  yet  it  predicts  the 
right  trends  and  magnitudes. 


IV  NATIVE  IMPERFECTIONS  GROWN  INTO  ALLOYS 


There  are  two  kinds  of  imperfections  we  shall  discuss:  vacancies 
and  dislocations.  At  present,  the  conclusions  are  tentative,  because  the 
detailed  theory  is  still  incomplete.  Vacancy  formations  energies  and  the 
melting  temperature  are  both  proportional  to  the  bond  energies  of  the 
constituents.  Thus,  one  would  expect  the  vacancy  density  just  below  the 
melting  point  to  be  about  the  same  for  all  materials  with  a  similar  coordi¬ 
nation  number;  in  our  semiconductor  case,  the  coordination  number  is 
four.  Dislocation  energies  are  mostly  owing  to  long-range  bond  angle  dis¬ 
tortions9  and  accordingly  are  sensitive  to  the  sheer  coefficients,  which 
vary  roughly  as  for  many  compounds.1  Based  on  this  general  argu¬ 
ment,  one  would  expect  to  find  fewer  dislocations  in  materials  with  shor" 


bond  lengths 
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Dislocations  appear  in  bulk  grown  materials  through  a  sequence  of 
steps.10  We  shall  only  discuss  the  simplest  case,  where  there  are  no 
mechanical  stresses  on  the  growing  material.  In  the  temperature  gra¬ 
dient  behind  the  growth  front,  there  is  usually  a  supersaturated  vacancy 
content.  These  excess  vacancies  may  diffuse  to  surfaces  where  they  do 
no  harm,  or  they  may  condense  into  vacancy  clusters  that  can  then 
organize  into  dislocation  loops.  If  a  loop  is  not  parallel  to  the  growth 
front,  it  can  climb  toward  the  front  by  absorbing  more  vacancies.  If  it 
climbs  faster  than  the  growth  progresses,  then  the  loop  will  reach  the 
front.  It  can  then  spread  and  grow  along  with  the  crystal.  Observed  from 
the  growth  surface,  such  a  loop  will  appear  to  be  two  edge  dislocations.  If 
this  dislocation  formation  process  can  be  inhibited  in  any  of  its  steps,  the 
resulting  material  will  have  fewer  dislocations. 

Zn  in  Zni-jCdjTe  or  Zni_,HgxTe  will  serve  this  function  in  two  ways 
The  first  is  that  Zn  stabilizes  the  weaker  bond  in  these  alloys;  accord¬ 
ingly,  vacancy  formation  energies  in  the  alloy  should  be  larger  than  those 
in  the  pure  constituent  CdTe  or  HgTe.  The  melting  temperature  will  also 
be  raised  somewhat,  which  effect  "wins”  is  a  matter  still  to  be  decided 
The  second  effect  is  clear  cut:  The  shorter  bond  lengths  of  the  ZnTe  will 
introduce  stiff  struts  into  the  system  that  will  increase  the  dislocation 
energies  and  thus  should  inhibit  dislocation  loop  formation  and  subse¬ 
quent  climb  The  large  dislocation  density  reduction  found  by  Bell  and 
Sen  in  Zn,_xCdxTe  with  only  4  percent  Zn  supports  this  speculation.11 


V  CONCLUSIONS 


We  have  argued  that,  in  an  alloy,  the  bond  energies  and  bond  lengths 
are  modified  from  their  pure  crystal  values.  Bond  energies  can  be 
strengthened  or  weakened  in  an  alloy  depending  on  the  relative  sp3  -state 
energies  and  the  consequent  net  electron  shift  between  the  constituents. 
For  most  compounds,  the  average  bond  length  modifications  for  each 
type  of  bond  are  generally  small  compared  to  those  suggested  by 
Vegard’s  rule,  in  agreement  with  experiment.  However,  HgCdTe  is  an 
exceptional  case,  because  the  bond  lengths  for  HgTe  and  CdTe  are  nearly 
the  same  (by  accident).  The  bond  length  and  energy  shifts  also  depend 
sensitively  on  the  local  configuration  and  the  bonds  are  locally  strained. 
These  effects  must  be  taken  into  consideration  in  a  proper  theory  of  alloy 
vacancy  formation  energies,  phonon  frequencies,  dislocation  energies, 
etc. 


The  conclusions  most  relevant  to  the  initial  question  posed  in  the 
introduction  are  that  for  the  narrow  gap  alloys: 

•  Cd  destabilizes  the  already  weak  HgTe  bond  and  has  little  effect 
on  the  dislocation  energy. 

•  Zn  both  stabilizes  the  HgTe  bond  and  increases  the  dislocation 
energy 

•  Zn  also  decreases  dislocation  density  of  CdTe  to  improve  sub¬ 
strates 

•  HgZnTe  electron  and  hole  mobilities  are  about  the  same  as  those 
of  comparable  bandgap  HgCdTe 

•  HgZnTe  is  likely  to  be  a  superior  materia!  to  HgCdTe  for  detector- 
operating  ,n  the  8-12  jum  regime 
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FIGURE  CAPTIONS 


Fig.  1  Schematic  View  of  Various  Bond  Configurations. 

(a)  A(3)B(1)  configuration 

(b)  An  A(2)B(2)  configuration 

(c)  An  A(3)B(1)  configuration  embedded  in  a  virtual  crystal. 

Fig.  2  Bond  lengths  as  a  function  of  concentration  x  for  the  Hg!_xCdxTe 
alloy  in  different  configurations. 

The  short  dashed  lines  are  for  the  CdTe  bonds  and  the  long-short 
dashed  lines  represent  the  behavior  of  the  HgTe  bonds. 

Fig  3  Bond  energy  as  a  function  of  concentration  x  for  the  Hgj_xCdxTe 
alloy  in  different  configurations. 

The  short  dashed  lines  are  for  the  CdTe  bonds  and  the  long-short 
dashed  lines  represent  the  behavior  of  the  HgTe  bonds.  The  solid 
lines  are  the  concentration  weighted  averages  of  the  respective 
bond  energies 

Fig  4  Schematic  picture  of  a  potential  configuration  diagram  for  the 
HgTe  and  CdTe  bonds  in  an  Hg(3)Cd(l)  configuration. 

For  the  HgTe  bond  the  energy  and  position  of  the  bond  center  are 
calculated  by  Method  1  of  Sec  II- B  and  the  final  stretched  position 
and  ener£»  cv  Method  2  The  results  quoted  for  the  CdTe  bond  are 
those  deh — ' , n e d  by  Method  2 


Fig.  5  AC  bond  energy  shifts  for  the  n-VT  compounds  in  A(l)B(3) 
configurations  expressed  as  a  percentage  of  the  pure  crystal  bond 
energy. 

The  B  atom  is  designated  in  the  line  that  shows  the  shift.  The  pure 
crystal  AC  bond  energies  and  lengths  are  given  at  the  top  of  the 
figure.  These  numbers  were  calculated  by  Method  1  in  Sec.  II-B. 

Fig.  6  Theoretical  versus  experimental  hardness. 

The  two  values  at  the  top  and  bottom  of  the  vertical  lines  desig¬ 
nated  for  each  compound  were  calculated  for  two  different  possi¬ 
ble  dislocation  arrangements.  The  lower  one  corresponds  to  the 
minimum  energy  configuration  and  the  upper  one  to  a  likely 
higher  energy  configuration.  The  experimental  numbers  were 
taken  from  Ref.  2. 


Table  I 


Correlated  term  values  and  bond  energies 
for  three  II-VT  compounds. 

The  experimental  bond  lengths,  polarity, 
and  metallic  contribution  to  the  bond  energy 
are  also  plotted. 


Correlated  Term  Values  [eV] 


Atom 

s 

P 

sp3 

Zn 

-10.224 

-4.920 

-6.246 

Cd 

-9.611 

-4.784 

-5.991 

Hg 

-10.946 

-4.872 

-6.391 

Te 

-19.620 

-9.824 

-12.273 

Bond  Energies 


Compound  d  [A]* 


^metallic  [eV^  Eb  [eV] 


Fig  1 


Schematic  View  of  Various  Bond  Configurations. 

(a)  A(3)B(  1)  configuration 

(b)  An  A(2)B(2)  configuration 

(c)  An  A(3)B(  1)  configuration  embedded  in  a  virtual  crystal. 
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orientatTan,  probably  lies  somewhere  between  Hm,n  and  Hj.  For  an 
indenter  with  t?  =  tt/4,  we  have 

Hmin  =  0.0969  (11) 

and  Hi/Hmin=  2.39.  Harrison9  has  shown  that  one  contribution  to  the 
shear  coefficient  (actually  Cn  —  Cjg)  is 

G  =  2.38  joules/ m3  .  (12) 

m  d® 

where  fir  is  Planck’s  constant  over  2rr,  m  is  the  free  electron  mass 
tl^/m  =  1.220  xlO-38  [joules-m2]),  ac  is  the  covalence, 
<*c  =  Va/(Vi,  +  Vfa)1/2.  and  d  is  the  bond  length.  We  will  approximate  G 
by  this  expresssion;  however,  we  also  have  contributions  from 
C**  (>  C„  —  Cj2)  for  some  crystal  orientations.  C44  should  have  a  similar 
functional  dependence  on  d,  but  a  larger  multiplicative  constant  factor  3 
Combining  Eqs.  (10)  and  (12),  using  v  -  0.2,  and  changing  the  dimensions 
to  those  in  terms  of  which  experimental  hardness  numbers  are  cus¬ 
tomarily  quoted  gives 


H 


mm 


2.38  x  10 


kgm/  mm2 


(13) 


where  d  is  in  3  Values  calculated  from  Eq.  (13)  and  Hj  are  plotted 
against  experimental  results  in  Fig  2  for  a  number  of  semiconductors 


for  the  minimum-energy  configuration  and  a  tetrahedrally  bonded  semi¬ 
conductor.  Finally,  in  the  indicated  configuration, 


for  all  i,  and  j. 


(8) 


Inserting  these  expressions,  Eqs.  (5),  (6),  (7),  and  (8)  into  the  sums 
in  Eq.  (3)  and  retaining  only  terms  of  order  N*  yields 


One  can  also  get  a  number  for  the  hardness  of  a  dislocation 
configuration  in  which  all  the  material  is  pushed  along  the  same  glide 
plane,  e.g.  to  the  side,  to  the  bottom,  or  normal  to  the  face8  of  the 
indenter  (a  possibility  not  depicted  in  Fig.  1).  In  this  case,  the  factor  of  2 
in  front  of  Eq.  (3)  is  removed,  N  =  Nj,,  and  cp  -  n/2  --0.  Then  a  higher 
nonequilibrium  hardness  in  the  context  of  this  model  (denoted  Hj,  is 
obtained  (remember  that  some  terms  have  been  neglected): 


H, 


G  cot  t?  cot  '3 

3rr  ( 1  -i/)  [  V2 


(10) 


The  proper  answer  for  most  materials,  and  depending  on  crystal 


-  '-“-tv  We  now  encounter  our  first  surprise.  As  we  can  see  from  Eqs.  (4)  and 
'  (5),  the  energies  per  unit  length  Ey  and  Ej  have  comparable  magnitudes. 
Because  there  are  approximately  N2  terms  in  the  interaction  energy  sum, 
only  N  terms  in  the  formation  energy  sum,  and  N  is  a  large  number 
(N  y>  1),  the  interaction  energy  completely  dominates  the  hardness.  In 
fact,  N  is  typically  of  order  10®.  Hence,  terms  owing  to  screw  disloca¬ 
tions,  core  energies,  heat  dissipation  as  the  dislocations  propagate  to 
their  places,  and  other  effects  associated  with  the  first  term  are  unim¬ 
portant.  However,  care  must  be  taken  with  the  interaction  terms:  Eq. 
(3)  neglects  a  number  of  secondary  interactions,  some  positive  and  oth¬ 
ers  negative;  these  will  be  added  later.  The  principal  neglected  terms  are 
the  interactions  between  the  dislocations  in  the  different  regions  on  each 
side  (positive),  and  the  interactions  between  adjacent  side  and  bottom 
regions  (positive),  and  the  interaction  between  the  opposite  sides  in  the 
bottom  regions  (negative).  Comparison  of  the  results  with  experiment 
will  indicate  how  important  these  neglected  terms  are  likely  to  be.  The 
length  Lj  is: 


1 

•? 

4 


Li  = 


W 


N-i 

N 


i  =  1,2,3 . N 


(e) 


V 


This  is  the  largest  length  of  the  side-inserted  planes  and  its  choice  par¬ 
tially  accounts  for  interactions  between  the  otherwise  neglected  screw 
dislocations.  The  distance  r^  =  r*  —  rj  is  given  by 

r jj  =  2^2b  !  i  -  j  !  ;  i,j  =  1,2.3, ...,N  (7; 
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the  shear  coefficient  is  G,  the  Burger’s  vector  b  =  d/3  for  an  indentation 
along  a  <100>  axis,  the  range  of  the  elastic  deformation  of  a  dislocation  R 
is  taken  equal  to  W  (for  want  of  a  better  approximation),  r0  ~  d  is  the 
dislocation  core  radius,  the  Poisson  ratio  is  v  ~  0.2  for  most  semiconduc¬ 
tors,  and  Lj  is  the  length  of  the  ith  dislocation.  In  the  second  term.  Ey  is 
the  interaction  energy  per  unit  length  between  dislocations  i  and  j, 
assuming  they  have  parallel  glide  planes  and  their  Burger’s  vectors  have 
the  same  sign.  (This  neglects  terms  that  will  be  discussed  presently.) 
The  expression  for  Ey  is8: 


where  ry  is  the  separation  between  dislocation  i  and  j,  and  $?y  is  the  angle 
that  a  line  perpendicular  to  and  joining  the  dislocations  makes  with  the 
glide  plane,  as  shown  in  Fig.  1.  Because  the  various  dislocations  in  a 
region  have  different  lengths,  the  net  interaction  energy  is  approximated 
by  multiplying  the  energy  per  unit  length  by  the  length  of  the  shorter 
one.  The  one-half  multiplying  this  term  prevents  double  counting  and  the 
prime  on  the  sum  indicates  that  terms  with  i  =  j  are  to  be  omitted.  The 
upper  limits  on  the  sums  N  are  the  number  of  dislocations  in  one  region 
(side  or  bottom)  associated  with  one  edge.  For  the  minimum-energy 
configuration,  N  =  N.,/2.  The  four  that  multiplies  the  bracket  accounts 
for  the  four  sides,  and  the  two  for  the  two  regions  for  each  side. 
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minimum-energy  configuration.  For  the  arrangement  depicted  in  Fig.  1, 
it  always  costs  energy  to  position  a  second  dislocation  on  a  parallel  glide 
plane  to  one  already  present.  However,  the  magnitude  of  this  extra 
energy  can  be  minimized  and,  for  proper  configurations  of  the  disloca¬ 
tions,  there  are  attractive  forces  along  the  glide  planes  that  will  tend  to 
position  the  dislocations  into  the  minimum-energy  configuration.  The 
minimum  configuration  arises  when  half  the  atoms  from  the  indented 
volume  go  respectively  to  the  side  and  below  the  indenter.  Then,  in  both 
regions,  the  maximum  angle  made  between  successive  close-spaced  dislo¬ 
cation  lines  and  their  glide  planes  is  <p  =  n/2  —$/  2,  as  shown  in  the 
figure.  This  is  the  minimum  realistic  energy  configuration.  As  we.  shall 
demonstrate  shortly,  the  minimum-energy  configuration  is  evidently 
established  when  soft  materials  are  indented;  however,  hard  ones  are 
more  apt  to  align  in  a  configuration  closer  to  a  higher  energy 
configuration  where  all  the  atoms  go  either  to  the  sides  or  below  the 
indenter.  If  the  dislocations  are  separated  more  than  shown  in  Fig.  1, 
then  there  is  more  volume  of  strained  material  and  the  interaction 
energy  would  be  larger  still. 

An  approximate  expression  for  the  energy  required  to  indent  the 
material  is8: 

tx=2  4  Ey  |  minOi.Lj)  J  (3) 

i*l  c  i,j*l 

In  the  first  term,  using  an  isotropic  medium  approximation,  the  energy 
per  unit  length  to  form  an  edge  dislocation  is8 


accommodate  an  indentation  to  depth  h  is 


N„=  f  =  .  (2) 

where  t?  is  the  angle  between  the  normal  to  the  tip  of  the  indenter  and  a 
side.  Figure  1  also  shows  a  model  of  one  possible  configuration  of  the 
dislocations.  The  edges  of  the  extra  atom  planes  that  are  driven  from  the 
indented  volume  into  the  bulk  of  the  semiconductor  are  shown  as  lines 
terminated  by  dots.  The  dotted  ends  of  these  lines  are  the  positions  of 
the  dislocations,  which  are  perpendicular  to  the  plane  of  the  figure.  The 
planes  driven  to  the  sides  each  have  a  finite  extent  and  a  trapezoidal 
shape.  The  planes  driven  down  under  the  indenter  have  a  square  shape. 

We  shall  argue  that  much  of  the  work  done  on  the  indenter  goes  into 
the  energy  to  form  the  indicated  dislocation  configuration,  although  some 
certainly  goes  into  heat.  There  are  two  major  contributions  to  this  for¬ 
mation  energy.  The  first  is  the  energy  needed  to  generate  each  disloca¬ 
tion  as  an  isolated  entity,  and  the  second  is  the  interaction  energy  among 
these  dislocations.8  Because  the  interaction  term  dominates  H,  approxi¬ 
mations  made  to  simplify  the  first  term  are  relatively  unimportant.  The 
extra  planes  driven  to  the  sides  of  the  indentation  have  a  finite  extent, 
accordingly  (in  this  idealized  picture),  there  are  both  edge  dislocations  at 
their  base  and  screw  dislocations  associated  with  their  termination.  The 
square  planes  driven  below  the  indentation  have  edge  dislocations  around 
the  sides  and  screw  dislocations  at  the  corners  to  make  the  turns.  More¬ 
over,  there  are  interactions  among  the  dislocations,  which  can  produce  ^ 
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arising  from  an  integral  over  the  energy  of  small  bond  angle  distortions 
distributed  over  many  sites.0  The  elastic  contribution  is  typically  several 
times  larger  than  the  core  energy;  therefore,  we  neglect  the  core  terms 
in  this  treatment,  where  the  major  trends  are  being  sought.  The  disloca¬ 
tion  energy  per  unit  length,  Eq  is  proportional  to  the  product  of  a  shear 
coefficient  and  the  square  of  the  Burger’s  vector.8  The  magnitude  of  the 
Burger’s  vector  is  proportional  to  the  lattice  constant8;  hence, 
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Hardness  is  determined  by  applying  a  known  force,  F,  to  a  probe  of  a 
prescribed  shape  driving  it  into  the  surface  of  the  sample.8  The  area.  A, 
of  the  resulting  indentation  is  measured,  and  the  hardness  is  the  force 
per  unit  indented  area.  Many  dislocations  must  be  formed  to  allow  the 
probe  to  indent  the  semiconductor.  If  the  indenter  is  a  rectangular 
pyramid,  then  the  hardness  is 


tt  _  F  _  Fh  _  CT 
A  Ah  Ah 


(1) 


where  e 7  is  the  work  required  to  cause  the  indenter  to  penetrate  to  a 
depth  h.  A  side  view  of  the  indentation  in  a  cut  through  its  center  is  illus¬ 
trated  schematically  in  Fig  1  The  top  of  the  indentation  has  side  length 
W;  thus,  A  =  W2  The  Burger's  vector  has  magnitude  b,  proportional  to 
the  bond  length  d  The  number  of  dislocations  N;n  required  to 
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metallic  contribution  to  the  bond  energy  in  semiconductors  is  also  bond- 

angle-dependent3. 

The  shear  coefficients  (combinations  of  C44  and  Cu  —  in  the 
Schoenflies  notation)®  depend  on  crystal  orientation  and  (in  Harrison’s 
notation®)  are  proportional  to  V|/d5(V|  +  V|  )3/’2,  where  V2  *  d-2  is  the 
covalent  and  V3  is  the  ionic  energy.  The  metallic  interaction  modifies  the 
functional  dependence  of  the  shear  coefficient  on  V3  and  V3,  but  intro¬ 
duces  no  explicit  dependence  on  the  hopping  integrals,  denoted  Vj  by 
Harrison.3,9  The  actual  independent  variables  that  determine  the  pro¬ 
perties  of  solids  are  the  valence  electron  atomic  energy  levels  and  wave 
functions.  Although  the  wave  functions  are  complicated,  their  essence 
for  trying  many  bonding  properties  is  captured  in  the  bond  lengths. 
Thus  the  atomic  energies  and  bond  lengths  can  be  used  as  the  effective 
independent  variables  of  the  problem.3  The  ionic  energy  V3  and  the  hop¬ 
ping  integrals  Vj  are  functions  only  of  the  atomic  energies,  while  the 
covalent  energy  V2  is  a  function  only  of  the  bond  length.  Thus,  in  a  case 
where  the  covalent  energy  dominates  (V2  »  V3)  the  shear  coefficients 
are  proportional  to  d-5;  in  the  other  limit,  they  are  smaller  and  propor¬ 
tional  to  d-u.  This  sharp  dependence  on  bond  length  results  from  the 
fact  that  the  shear  coefficients  are  roughly  the  covalent  bond  energy  con¬ 
tribution  per  urn',  bond  volume.  In  a  pure  covalent  material,  the  bond 
energy  is  proportional  to  V2  (or  d-2),  and  the  bond  volume  is  «  d3;  hence, 
in  this  case,  the  shear  coefficient  varies  as  d-® 

There  are  two  contributions  to  the  dislocation  energy  per  unit 
length:  a  short-ranged  core  energy  and  an  elastic  long-range  energy 
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Dislocations  in  semiconductors  are  detrimental  to  device  function; 
they  serve  as  channels  for  impurity  migration  and  trapping,  which  causes 
nonuniform  doping  and  degrades  p-n  junctions.1  They  also  decrease  the 
material’s  resistance  to  plastic  deformation.  The  aim  of  this  paper  is  to 
provide  insights  into  the  underlying  physical  mechanisms  controlling 
dislocations  and  semiconductor  hardness,  and  then  to  suggest  strategies 
for  decreasing  dislocation  densities.  It  is  well-established  that  the  hard¬ 
ness  of  tetrahedrally  coordinated  semiconductor  materials -Group  IV,  in- 
V,  and  II-VI  compounds— exhibits  a  sharp  variation  with  their  near- 
neighbor  distance  d,  approximately  proportional  to  d“®  for  one  group  of 
seven  compounds.2  (The  lattice  constant,  a,  and  d  are  related  by 
d  =  aV3/4.)  Thus,  semiconductors  with  small  lattice  constants  tend  to 
be  harder  materials.  These  same  materials  have  larger  stiffness 
coefficients3  and  have  fewer  dislocations  in  as-grown  crystals.4-6  We  now 
interrelate  these  observations  and  some  other  more  recent  measure¬ 
ments,  through  an  extension  of  Harrison’s  theory  of  the  bonding  of  these 
semiconductors  3 

The  major  reason  for  the  properties  just  described  can  be  under¬ 
stood  qualitatively  as  follows.  These  materials  have  a  low  coordination 
number  with  only  four  near  neighbors.  There  are  three  contributions  to 
their  bond  energies.  The  most  important  one  is  the  covalent  energy  that 
results  from  the  interactions  between  the  sp3  hybrids  pointing  toward 
one  another  on  adjacent  sites.  The  strength  of  this  interaction  is  direc¬ 
tionally  dependent  and  energy  is  required  to  cause  a  bond  angle  distor¬ 
tion.  Secondly,  the  direct  ionic  contribution  to  the  bond  tension  is  radial 
(i.e.  nondirectional)  and  adds  little  to  the  shear  coefficients  Thirdly,  the 
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DISLOCATION  ENERGIES  AND  HARDNESS  OF  SEMICONDUCTORS 

A.  Shef ,  SRI  International,  Menlo  Park,  CA  94025 
A.-B.  Chen,  Auburn  University,  Auburn,  AL  36849 
W.E.  Spicer,  Stanford  University,  Stanford,  CA  95305 

ABSTRACT 

* 

The  dislocation  energies  and  hardness  of  semiconductors  are  calcu¬ 
lated  by  an  extension  of  Harrison's  method.  It  is  demonstrated  in  agree¬ 
ment  with  experiment  that  dislocation  energies  per  unit  length  are  pro¬ 
portional  to  d-3  to  d-9.  where  d  is  the  bond  length,  and  hardness  is  pro¬ 
portional  to  d-5  to  d-n.  The  hardness  is  related  to  the  interaction  ener¬ 
gies  among  dislocations.  It  is  argued  that  dislocation  densities  of  semi¬ 
conductors  will  be  reduced  if  they  are  alloyed  with  a  second  constituent 
that  has  a  shorter  bond  length.  Experimental  evidence  supporting  this 
strategy  is  noted. 
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Fig.  6  Theoretical  versus  experimental  hardness.  / 

-  3 

The  two  values  at  the  top  and  bottom  of  the  vertical  lines  desig-  • 

■ 

nated  for  each  compound  were  calculated  for  two  different  possi-  .• 

ble  dislocation  arrangements.  The  lower  one  corresponds  to  the 

minimum  energy  configuration  and  the  upper  one  to  a  likely  w 

- 

higher  energy  configuration.  The  experimental  numbers  were  ; 

taken  from  Ref.  2. 
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Fig.  5  AC  bond  energy  shifts  for  the  II-V]  compounds  in  A(l)B(3) 
configurations  expressed  as  a  percentage  of  the  pure  crystal  bond 
energy 

The  B  atom  is  designated  in  the  line  that  shows  the  shift.  T^e  pure 
crystal  AC  bond  energies  and  lengths  are  given  at  the  top  of  the 
figure.  These  numbers  were  calculated  by  Method  1  m  Sec.  tl-B 
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Fig.  4  Schematic  picture  of  a  potential  configuration  diagram  for  tne 
HgTe  and  CdTe  bonds  in  an  Hg(3)Cd(l)  configuration. 

For  the  HgTe  bond  the  energy  and  position  of  the  bond  center  are 
calculated  by  Method  1  of  Sec.  II-B  and  the  final  strecthal  position 
and  energy  by  Method  2.  The  results  quoted  for  the  CdTe  bond  are 
those  determined  by  Method  2 
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CONFIGURATION 
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Fig.  3  Bond  energy  as  a  function  of  concentration  x  for  the  Hgj_xCdxTe 
alloy  in  different  configurations. 

The  short  dashed  lines  are  for  the  CdTe  bonds  and  the  long-shor 
dashed  lir.es  represent  the  behavior  of  the  HgTe  bonds.  The  sohc 
lines  are  the  concentration  weighted  averages  of  the  respeciivt 
bond  energies 


Figured  has  the  theoretical  Ham  and  Hj  values  connected  by  arrows 
from  ffmin  to  Hj  for  each  compound,  plotted  as  a  function  of  the 
corresponding  experimental  values.8,10  If  the  theory  were  perfect  and 
the  experimental  values  were  accurate,  the  points  would  fall  on  the  indi¬ 
cated  unity  slope  line.  Several  conclusions  can  be  drawn.  First,  the 
order  of  magnitude  of  the  predicted  and  measured  values  are  the  same,  a 
result  obtained  with  no  adjustable  parameters  in  the  theory.  Secondly, 
the  trends  from  one  compound  to  another  are  properly  given  by  the 
theory.  Although  the  values  are  generally  too  small,  they  fit  the  soft 
materials  better,  and  the  Hi  values  fit  the  harder  materials  better.  As  we 
have  pointed  out,  some  terms  have  been  neglected  that  will  very  likely 
increase  H^n  (and  Hi  ),  but  they  should  be  smaller  than  those  retained. 
Thirdly,  from  Eqs.  (1)  and  (3),  H  is  given  in  a  rough  but  revealing  approxi¬ 
mation  by  4N2EyL/  Ah,  where  the  average  dislocation  length  L  =  W/2  and 
Eij  is  the  average  dislocation  pair  interaction  energy  per  unit  length. 
Notice  that  N  (or  W)  cancels  from  this  expre  :  .on;  thus,  H  is  independent 
of  W  (or  F)  and  therefore,  H  is  truly  a  measure  of  the  properties  of  the 
material.  This  result  would  not  be  found  if  the  dislocation  energies  («  N) 
were  to  dominate  H  rather  than  the  pair  interaction  energy  (°<  N2). 
Finally,  the  Berger’s  vector  cancels  from  the  leading  term  and  appears 
only  in  the  argument  of  the  logarithm  in  Eq.  (5).  Thus,  the  answers  are 
also  insensitive  to  its  choice. 

Dislocations  are  often  found  in  materials  as  they  are  grown.  Their 
density  is  determined  by  the  thermal  and  mechanical  stresses  to  which 
they  are  srbjected  in  the  growth  process.  A  dislocation  constitutes  a 
metastable  excitation  relative  to  the  perfect  crystal  ground  state.  At  the 


elevated .  growth  temperatures  and  temperature  gradient  behind  the 
growth  front,  the  number  of  dislocations  present  is  controlled  by  the 
relative  rate  at  which  vacancies  anneal  or  condense  into  dislocations.11 
The  dislocation  formation  rate  will  be  slower  in  a  material  grown  at  the 
same  temperature  if  Eq  is  higher.  If  an  alloy  is  formed  from  a  material  of 
interest  and  a  second  constituent  with  a  shorter  bond  length,  one  expects 
the  average  bond  energy  (and  thus  both  the  melt  temperature  and 
vacancy  formation  energy)  to  increase  proportional  to  a  low  inverse 
power  of  the  average  bond  length.9  Hence,  the  equilibrium  vacancy  den¬ 
sity  just  below  the  melting  point  tends  to  be  the  same  in  lowest  order  for 
all  materials,  independent  of  the  bond  lengths  of  the  constituents.  How¬ 
ever,  the  shear  coefficient  and  dislocation  energy  per  unit  length  will 
increase  with  much  higher  inverse  powers  of  the  bond  length.  Conse¬ 
quently,  dislocation  densities  should  be  reduced  in  such  alloys  relative  to 
those  found  in  the  longer  bond  length  pure  constituent.  This  expectation 
is  confirmed  in  the  recent  work  reported  on  Zn^Cd^Te  bulk  material.4 
The  best  CdTe  that  has  ever  been  grown  has  dislocation  densities  in 
excess  of  5  x  105  cm-2.  The  addition  of  only  4%  Zn  reduced  the  disloca¬ 
tion  count  to  less  than  5  x  104  cm-2.  The  ZnTe  bond  length  is  2.643  A, 
while  the  CdTe  bond  length  is  2.905  A,  a  6%  difference.  This  6%  difference 
in  bond  length  translates  into  a  2%  difference  in  the  dislocation  energy 
per  unit  length  for  1-x  =  0.04.  Dislocation  energies  per  unit  length  are 
typically  10  eV  per  lattice  spacing;  accordingly,  a  2%  increase  can  be 
expected  to  slow  their  formation  rate  considerably. 


:  ~The?lrgument  just  presented  naturally  leads  to  a  strategy  for 
'  decreasing dislocations  in  other  semiconductors.  If  an  alloy  is  made  of 
the  material  of  interest  with  another  compound  with  a  shorter  bond 
length,  then  the  dislocation  density  should  be  reduced.  For  example,  this 
suggests  that  the  addition  of  a  small  amount  of  GaP  (d  =  2.359  A)  may 
significantly  reduce  the  dislocation  density  of  bulk  grown  GaAs  (d  =  2.448 
A).  It  has  been  demonstrated  that  the  addition  of  approximately  1%  GaN 
(d=  1.948  A)5,8  or  of  a  1018  cm-3  BAs  (d  =  2.069  A)  concentration12  to 
GaAs  can  yield  a  large  volume  of  dislocation-free  material.  An  InAs  addi¬ 
tive  with  its  longer  bond  length  (d  =  2.623  A)  serves  the  same  function 
indirectly,  by  causing  GaAs  bonds  in  its  neighborhood  to  be  compressed. 
This  indirect  mechanism  should  be  less  effective  than  substituting  short- 
bond-length  additives. 

We  have  demonstrated  that  the  dislocation  energies  and  hardness  of 
tetrahedrally  bonded  semiconductors  are  rapid  functions  of  the  recipro¬ 
cal  of  the  bond  length.  This  rapid  d  dependence  of  dislocation  energies 
provides  a  rationale  for  the  dramatic  decrease  of  the  dislocation  density 
in  bulk  grown  Zno.04Cdo.98Te  material  relative  to  that  found  in  CdTe,  and 
suggest  means  for  accomplishing  the  same  ends  in  other  materials. 
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FIGtJRE  CAPTIONS 


Fig.  1  A  schematic  representation  of  an  ideali2ed  minimimum-energy 
dislocation  configuration  produced  by  a  square-cross-section 
indenter  in  a  hardness  measurement. 

W  is  the  side  length  and  h  is  the  depth  of  the  indentation.  In  this 
ideal  case,  half  the  material  from  the  indented  region  is  displaced 
along  the  glide  planes  (indicated  by  the  dashed  arrows)  to  the 
sides  and  half  is  displaced  below  the  indenter. 

Fig.  2  The  theoretical  vs.  experimental2,10  hardness  of  several  semicon¬ 
ductors. 

The  two  theoretical  values  for  each  semiconductor  are  connected 
by  an  arrow  from  to  H*.  Perfect  agreement  would 

correspond  to  points  being  on  the  unity-slope  solid  line. 
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Fig.  1  A  schematic  representation  of  an  idealized  minimimum-energy 
dislocation  configuration  produced  by  a  square-cross-section 
indenter  in  a  hardness  measurement. 

W  is  the  side  length  and  h  is  the  depth  of  the  indentation.  In  this 
ideal  case,  half  the  material  from  the  indented  region  is  displaced 
along  the  glide  planes  (indicated  by  the  dashed  arrows)  to  the 
sides  and  half  is  displaced  below  the  indenter. 
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H[kjm  |mm)2]  —  txp*nm*nt 


Fig.  2  The  theoretical  vs.  experimental2,10  hardness  of  several  semicon¬ 
ductors. 

The  two  theoretical  values  for  each  semiconductor  are  connected 
by  an  arrow  from  Hmi,  to  Hi-  Perfect  agreement  would 
correspond  to  points  being  on  the  unity-slope  solid  line. 
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